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The  objective  of  this  research  is  to  develop  techniques  and  procedures  for 
computing  the  statistical  characteristics  of  some  electromagnetic  scattering 
phenomena  from  radar  targets.  Radar  measurements  are  perturbed  by  electromagnetic 
scattering  phenomena  induced  by  complex  targets  including  radar  cross  section 
scintillation  and  glint  or  position  error.  The  statistical  characteristics  of  these 
phenomena  may  be  used  to  evaluate  radar  system  performance,  develop  optimum  tracking 
systems,  and  develop  stochastic  models  to  use  in  system  simulations. 

The  procedures  developed  use  numerical  integration  to  compute  the  moments  of 
the  probability  density  functions  and  the  characteristic  function  method  to  compute 
the  covariance  functions.  As  part  of  the  research,  a  deterministic  model  of  the 
electromagnetic  scattering  from  M  bodies  has  been  developed  and  applied  to  an 
ellipsoidal  model  of  a  common  target  drone.  This  model  is  used  as  the  electromag* 
netic  scattering  input  to  the  analytic  computation  of  the  statistical  characteris¬ 
tics  and  to  the  simulation  which  is  used  to  check  the  general  accuracy  of  the 
analytic  solutions. 

The  statistical  characteristics  are  computed  at  three  points  on  a  typical 
flight  path  of  interest  for  three  different  sets  of  target  flight  characteristics. 
These  statistical  characteristics  are  conditional  on  the  flight  path,  i.e.,  they 
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are  functions  of  the  flight  path.  The  agreement 
between  the  analytic  computation?)  and  simulations 
is  excellent.  The  sensitivity  of  the  electromag¬ 
netic  scattering  statistics  to  the  flight  motion 
statistics  is  generally  as  expected.  The  data 
from  the  simulations  are  very  similar  in  nature 
to  dynamic  measurements. 
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1 .  INTRODUCTION 

1.1  Introduction  to  the  Radar  Problem 

Many  modern  radar  systems  are  capable  of  tracking  targets  In  four 
dimensi'.,s:  range,  azimuth  angle,  elevation  angle,  and  doppler 
frequency.  There  is,  however,  a  point  at  which  target  induced  effects 
become  significant  contributors  to  system  measurement  errors.  Angle 
tracking  is  perturbed  by  an  effect  called  glint;  range  tracking  and 
other  functions  are  perturbed  by  cross-section  scintillation;  and 
doppler  tracking  is  perturbed  by  apparent  target  phase  center  motion 
not  related  to  translational  motion. 

A  pulse  doppler  radar  system  transmits  a  pulse  of  electromagnetic 
energy  of  known  amplitude  and  phase  at  a  known  time.  The  signal 
reflected  from  a  target  is  received  by  the  radar,  and  is  processed  to 
give 

(1)  The  target  range 

(2)  The  target  size  or  radar  cross-section  (RCS) 

(3)  The  target  azimuth  angle 

(4)  The  target  elevation  angle 

(5)  The  target  phase  angle. 

The  errors  involved  in  these  measurements  can  be  divided  into  two 
basic  types.  Type  I  errors  are  those  caused  by  noise  and  imperfect 
signal  and  data  processing.  These  errors  exist  even  for  a  perfect 
target  such  as  a  sphere.  Type  II  errors  are  those  introduced  by  the 
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the  target.  A  complex  target  such  as  an  aircraft*  or  even  two  spheres, 
has  RCS  and  scattering  phase  which  are  functions  of  the  aspect  angles 
and  reflect  or  scatter  a  wavefront  which  is  distorted  or  nonspherical. 

The  Type  I  errors  have  been  investigated  and  extensively  studied 
[1,  2,  3] .  The  Type  II  errors  have  had  less  development  than  the  Type 
I  errors.  The  primary  reasons  for  this  are  the  complexity  of  the 
problem  and  the  dependence  of  the^effects  on  the  particular  body  under 
consideration.  Unlike  thermal  noise  which  can  readily  be  modeled  as 
white  noise,  the  target  induced  effects  depend  on  the  target,  its 
aspect  angles  and  its  motion.  The  following  is  a  summary  of  the  errors 
of  interest  and  an  indication  of  the  relationship  of  the  two  types  of 
errors , 


The  equation  for  power  received  by  a  monostatic  radar  is  given  by 
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(1.1) 


where 


Pt  s  transmitted  power 
G  =>  antenna  gain 
\  a  radar  wavelength 
R  e  range  from  radar  to  target 
L  ■  numerical  factor  to  account  for  losses 


S  a  target  RCS. 

The  measurement  of  target  range  is  generally  limited  by  the  signal-to- 
noise  ratio  which  depends  on  received  power.  The  rms  range  error 
caused  by  thermal  noise  is  given  by  ( 3] : 


where 


c  ■  velocity  of  light 
t  =  transmitted  pulse  width 
B  «  receiver  bandwidth 
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E  =  received  energy 
Nq  a  receiver  noise  power  spectrum. 

The  target  induced  errors  are  manifested  as  changes  in  the  received 
energy,  E,  (sometimes  exceeding  50  decibels  peak  to  peak)  caused  by 
changes  in  apparent  target  RCS;  it  is  apparent  from  (1.2)  that  the  two 
errors  are  interrelated.  Similarly,  the  target  induced  errors  affect 
the  angular  position  measurements  and  the  target  phase  angle  (doppler 
frequency)  measurements . 

All  angle  tracking  systems  are  essentially  phase-front  measuring 
devices  and  have  been  shown  to  be  equivalent  to  each  other  in  perfor¬ 
mance  [4],  Only  amplitude  monopulse  systems  will  be  discussed  herein 
since  they  are  predominant  in  the  modern  generation  of  radar  systems. 

A  dual-plane  monopulse  radar  generally  has  four  squinted  beams  which 
are  combined  to  give  a  sum  channel  and  two  angle  error  channels  [2] . 

An  ideal  target  on  boresight  will  give  zero  outputs  from  the  angle 
tracking  channels.  The  target  induced  error  is  manifested  as  a  tilt 
in  phase  front  which  gives  an  apparent  target  error  signal.  The 
angular  error  signal  is  inversely  proportional  to  range.  The  effects 
of  glint,  then,  are  most  significant  at  close  ranges  such  as 
encountered  in  homing  missile  systems.  The  theoretical  rms  angular 
error,  caused  by  thermal  noise  for  a  uniformly  illuminated  antenna, 
is  given  by  [1] 
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0.628  e, 
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(1.3) 


where  8^  is  the  half-power  beam  width.  Again  the  effect  of  variation 
in  E  should  be  noted. 

Doppler  tracking,  in  a  pulse  doppler  radar,  is  accomplished  by 
measuring  the  phase  shift  of  the  RF  signal  from  transmission  to 
reception,  modulo  2fl  radians.  The  rms  phase  error  caused  by  thermal 
noise  is  given  by  [3] 


e 

<P 


(1.4) 


In  addition  to  this  error,  there  is  an  apparent  target  phase  center 
motion  which  is  aspect  angle  dependent.  Target  motion,  other  than 
purely  radial,  may  cause  a  random  modulation  of  the  target  phase  range. 

However,  as  outlined  previously,  the  major  interdependence  of  the 
two  types  of  errors  are  caused  by  the  change  in  RCS  and  so  in  the 
received  energy,  E.  Therefore,  it  is  possible  to  separate  the  two 
types  of  errors  in  a  convenient  manner.  It  will  be  assumed  hereafter 
that  the  only  interdependence  is  that  induced  by  the  RCS.  Thus, 
modeling  of  the  target  induced  effects  need  not  consider  the  radar 
characteristics,  except  as  delineated  below,  and  the  radar  receiver 
model  need  not  be  modified  to  consider  the  Type  II  errors. 

The  phenomena  of  amplitude  scintillation,  glint,  and  doppler 
scintillation  are  different  manifestations  of  the  vector  summation  of 
the  electromagnetic  waves  reflected  and  refracted  by  complex  targets. 
Although  basically  target  induced,  these  phenomena  are  affected  by 
some  of  the  characteristics  of  the  observing  radar  such  as  resolution, 
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polarization,  wavelength,  tracking  implementation,  and  type  (monostatic 
or  bistatic*) .  All  but  very  few  targets  have  RCS,  glint,  and  phase 
centers  which  are  sensitive  to  wavelength  and  polarization.  These 
characteristics  should  be  included  in  any  complete  model  of  the  target 
phenomena,  but  not  all  are  within  the  scope  of  this  work. 

Although  bistatic  radars  are  of  some  interest,  especially  in  semi¬ 
active  homing  systems,  the  discussion  and  proposed  research  will  be 
limited  to  monostatic  systems.  It  is  assumed  that  the  range  resolution 
is  such  that  the  entire  target  is  always  in  the  range  resolution  cell. 

Similarly,  it  will  be  assumed  that  the  target  is  smaller  than  th£ 

» 

< 

angular  resolution  and  is  always  in  the  angular  resolution  cell.  The 
range  of  radar  frequencies  will  be  limited  to  those  most  commonly  used 
for  search  and  tracking  systems,  i.e.,  L-band  through  X-band  GHz 
through  10  GHz) . 

It  would,  at  first,  appear  to  be  of  value  to  have  a  complete 
deterministic  model  of  the  particular  target  of  interest  to  the  radar 
system.  Further  consideration,  however,  indicates  that  this  is  not 
a  very  practical  goal.  The  target  appearing  in  the  radar  field  of  view 
and  the  target  aspect  angles  will  seldom  be  known,  except  perhaps  by 
the  ability  to  classify  the  target  in  general  terms  such  as  rotary 
wing;  fixed  wing,  propeller  driven;  fixed  wing,  jet  propelled;  etc. 

It  would  not  be  practical  to  store  complete  target  data  for  all  targets 
of  interest  in  the  radar  system  computer  even  if  it  was  possible  to 


A  monostatic  radar  has  its  transmitting  and  receiving  antennas 
collocated  or  common  whereas  a  bistatic  radar  may  have  its  trans¬ 
mitting  and  receiving  antennas  separated  by  large  distances. 
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generate  such  data.  The  most  obvious  and  most  practical  approach  is  to 
model  the  problem  statistically. 

The  objective  of  this  research  has  been  to  develop  techniques  to 
determine  important  statistical  properties  of  target  induced  scintilla¬ 
tion  and  glint. 

1.2  Review  of  Past  Work 

A  brief  review  of  the  history  and  evolution  of  radar  systems  and 
the  development  of  models  will  help  to  gather  perspective  on  the 
problem  and  the  status  of  work  in  the  area.  The  review  presented  is 
necessarily  brief  and  oriented  toward  the  problem  under  study. 

References  [1,  2,  3,  5,  6,  7,  and  8]  are  general  treatises  on  modern 
radar  systems. 

The  radar  systems  operating  at  the  start  of  World  War  II  were  used 
primarily  for  early  warning  and  ranging  since  they  had  poor  angular 
accuracy.  The  war  spurred  research  and  development  on  radar  systems, 
especially  radar  fire  control  systems  for  antiaircraft  systems.  The 
SCR  584,  introduced  early  in  the  war,  was  the  first  operational  fire 
control  system  developed  by  the  U.S.  which  required  no  assistance  from 
optical  angle  trackers.  It  could  track  targets  in  three  dimensions: 
azimuth  angle,  elevation  angle,  and  range.  Since  World  War  II  the 
development  and  evolution  of  radar  systems  have  continued  at  an  almost 
fantastic  pace.  The  types  of  radar  systems  have  proliferated  due  to 
differing  operating  requirements  and  constraints. 

Amplitude  scintillation  was  noted  in  World  War  II  when  target 
fading  was  frequently  observed.  It  was  not  until  angle  tracking  systems 
improved  significantly  that  glint  was  first  noticed.  At  first  it  was 
ascribed  as  an  effect  of  amplitude  scintillation  since  large  values  of 
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glint  occurred  only  with  low  signal  levels.  Later  it  was  determined  to 
be  a  separate  phenomenon.  The  relationships  between  glint  and 
amplitude  scintillation  had  to  wait  even  longer  before  it  was  clearly 
established.  Zt  is  now  accepted  that  the  three  phenomena  of  amplitude 
scintillation,  glint,  and  doppler  scintillation  are  different  manifes¬ 
tations  of  the  same  target  induced  effects. 

The  development  of  target  models,  both  analytical  and  statistical, 
parallels  the  evolutionary  growth  in  sophistication  and  performance  of 
radar  systems.  Amplitude  scintillation  and  its  effects  on  target 
detection  have  been  extensively  treated  from  a  statistical  point  of 
view.  Marcum  [9]  developed  the  theory  of  detection  probabilities  of 
constant  cross  section  targets  in  additive  white  noise  for  search  type 
radars,  considering  the  effects  of  pulse  integration,  scanning  effects, 
etc. 
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Swerling  [10,  11]  advanced  the  theory  by  considering  fluctuating 
targets.  The  four  cases  originally  considered,  which  encompassed  most 
targets,  are  as  follows: 

(1)  Case  1  -  The  target  is  assumed  to  be  constant  for  time-on- 
target  for  a  single  scan  but  fluctuates  randomly  on  a  scan-to-scan 
basis.  The  probability  density  function  of  the  target  cross  section, 
S,  is  assumed  to  be: 

fs<»)-  exp  ^  ,  s  &  0  (1.5) 

where 


s  »  E  (S) 
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(2)  Case  2  -  The  target  fluctuations  are  random  on  a  pulse-to- 
pulse  basis.  The  probability  density  function  is  the  same  as  for 
Case  1. 


(3)  Case  3  -  The  target  is  assumed  to  be  constant  for  the  time- 
on-target  of  a  single  scan  but  fluctuates  on  a  scan-to-scan  basis  with 
the  cross  section  having  a  probability  density  function  given  by 


exp  j  ,  s  ^  0  .  (1.6) 

(4)  Case  4  -  The  target  fluctuations  are  random  on  a  pulse- to- 
pulse  basis  with  the  same  probability  density  function  as  Case  3. 

These  distributions  are  chi-square  distributions  with  2  and  4 
degrees  of  freedom  for  Cases  1  and  2  and  Cases  3  and  4,  respectively. 

Swerling's  work  has  found  wide  acceptance  in  the  field,  especially 
since  the  results  were  presented  in  convenient  graphical  form.  Swerling 
[12]  has  investigated  additional  fluctuating  models  including  log- 
normally  distributed  targets.  Heidbreder  and  Mitchell  [13]  have  also 
investigated  detection  probabilities  for  log-normally  distributed  tar¬ 
gets.  A  log-normally  distributed  target  is  one  which  has  a  cross 
section  whose  logarithm  is  normally  distributed.  The  density  function 
is  given  by 


fs(s) 


as 


2  it 


exp 
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s  &  0 
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where  p  and  a  are  the  mean  and  variance  of  In  S,  respectively. 
Sponsler  [14]  investigated  the  track-while-scan  problem  encountered 
with  mechanically  scanned  search  radars.  Using  a  first-order,  two- 
state  Markov  process,  Sponsler  established  bounds  on  the  scan-to-scan 
correlation  coefficient  and  derived  the  Kolomogorov  differential 


equations  for  a  nonstationary,  continuous  parameter  Markov  process 
which  could  be  used  to  simulate  radar  data  for  arbitrary  flight  paths . 
The  equations  are  a  function  of  the  blip-to-scan  ratio,  i.e.,  time-on- 
target  to  scan  time  ratio. 

Weinstock  [15]  investigated  target  models  for  missile  and 
satellite  shapes  and  concluded  that  the  chi-square  models  give  poor 
representation  of  the  tails  of  the  distributions.  Weinstock  also  con¬ 
cluded  that  if  one  was  constrained  to  using  chi-square  models  the 
median  rather  than  the  mean  should  be  specified. 

The  U.S.  Naval  Research  Laboratory  has  been  investigating  glint 
since  1947.  The  efforts  have  been  devoted  primarily  to  understanding 
and  explaining  the  underlying  causes  of  the  phenomenon.  The  first 
published  work  on  statistical  modeling  of  glint  was  by  Delano  [  16 17]. 
Considering  mathematically  simple  arrays  of  point  scatters,  Delano 
computed  the  statistical  distribution  of  the  angular  errors.  For 
example,  considering  a  single  angle  tracking  channel  and  a  linear  array 
of  scatterers  of  statistically  independent  amplitude  and  phase,  all 
scatterers  having  approximately  the  same  mean  value,  Delano  derived  the 
result  that  the  apparent  target  center  is  outside  the  actual  target 
13.4  percent  of  the  time.  This  result  was  of  great  significance  and 
has  been  quoted  many  times,  often  out  of  context. 

Muchmore  [18,  19]  investigated  amplitude  scintillation  spectra 
using  Delano's  models.  Although  Muchmore's  and  Delano's  work  was 
criticized  by  Peters  and  Weimer  [20,  21,  22]  because  of  the  simplicity 
of  the  models,  the  spectra  obtained  were  reasonably  similar  to 
empirical  results.  Development  of  statistical  glint  models  has  not 
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progressed  significantly  since  these  two  papers,  probably  due  in 
part  to  the  extreme  complexity  of:  the  problem, 

Meade  et  al.  [,?3]  investigated  the  two  body  model,  often  called 
the  dumbbell  model,  tnd  derived  the  result  that  the  apparent  target 
position,  0T>  is  given  by 


»  /  v-  *2  \ 

1  2  \  1  +  k2  +  2k  cos  q>  / 


(1.8) 


where 


0q  a  angular  separation  of  the  targets 
k  «=  cross  section  ratio  of  the  targets 
c p  =  relative  phase  angle  of  the  response 
of  the  tarots. 

Meade's  results  indicate  that  the  ;ngle  error  can  go  to  infinity. 
Ostrovityanov  [ 24]  corrected  this  errot  by  noting  that  the  assumption 
tan  0  «  0  was  implicit  in  Meade's  derivat:on.  Making  the  correction, 
one  finds  that  0  is  bounded  by  ±jf/2  which  is  more  reasonable.  Howard 
[  25]  interpreted  glint  as  a  distortion  or  tii*‘  in  the  wavefront  and 
demonstrated  means  of  computing  glint  as  a  function  of  aspect  angle 
for  N-body  targets,  where  N  is  any  finite  number,  Gubonin  [26]  used 
Howard's  interpretation  of  glint  and  derived  glint  statistics,  arriving 
at  essentially  the  same  results  as  Delano.  Gubonin  die'-,  however, 
avoid  the  implied  approximation  of  tan  0  «  0  and  obtained  more  general 
results . 

Lindsay  [27]  expanded  on  the  phase  front  approach  and  concluded 
that  glint  and  amplitude  and  doppler  scintillation  can  be  exp*'.  i:aed  by 
considering  Vcp,  i.e.,  the  normal  to  the  phase  front,  and  its  relation 
to  the  antenna  beam  axis.  Dunn  and  Howard  [28]  concluded  basically 
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the  same  fact  Independently,  and  demonstrated  that  the  phenomena  can 
be  explained  by  considering  the  Poynting  vector  of  the  reflected  or 
scattered  wave. 

Modeling  of  large  complex  radar  targets  such  as  aircraft  is  at 
best  an  extremely  difficult  and  rather  inexact  task  [29],  and  the 
techniques  used  are  guided  to  seme  extent  by  the  specific  objective 
of  the  researcher.  Little  or  no  effort  has  been  devoted  to  target 
modeling  as  proposed  herein.  Previous  work  has  been  devoted  to 
investigating  or  modeling  one  or  more  of  the  phenomena  such  as  Borison 
130]  in  doppler  scintillation;  Swerling  [10,  11,  12],  Weinstock  [15], 
Heidbreder  and  Mitchell  [13],  Edrington  [31]  and  Sponsler  [14]  in 
amplitude  scintillation;  Muchmore  [18]  in  scintillation  spectra; 

Delano  [16],  Howard  [25],  Gubonin  [26],  Besalov  and  Ostrovityanov  [32], 
Di— <u  and  Howard  [29,  33]  ,  and  Sims  and  Graff  134,  35]  in  glint;  or  in 
the  analysis  of  the  relationships  among  the  phenomena  such  as  work  by 
I.indsay  [27]  and  Dunn  and  Howard  [28]  . 

RCS  modeling  is  the  most  developed  area  of  radar  target 
scattering.  The  IEEE  devoted  an  entire  issue  of  the  Proceedings  (Vol. 
53,  No.  8,  August  1965)  to  radar  reflectivity.  Included  in  the  papers 
was  a  rather  extensive  bibliography  [36]  of  work  in  the  field  through 
1964.  A  bibliography  of  more  recent  works  may  be  found  in  [37] . 
Techniques  in  RCS  modeling  range  from  the  exact  solution  for  spheres 
by  Mie  [38]  to  the  empirical  approximations  by  Crispin  and  Maffett 
[39].  Three  comprehensive  works  in  the  field  have  been  published  in 
recent  years  [40,  41,  42] .  Some  of  the  more  interesting  work  in  radar 
scattering  recently  include  the  works  of  Ross  [43,  44],  Oshiro,  et  al., 
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[45,  46,  47],  Moll  and  Seecamp  [48],  Ryan  <49,  50],  Keller  [51,  52, 

53],  Bechtel  [54],  Uslenghi  and  Lee  [55] ,  and  Pierson  and  Clay  [56] . 

Electromagnetic  (EM)  scattering  has  been  categorized  into  three 
regions:  Rayleigh,  resonance,  and  optical,  depending  on  the  character¬ 
istic  dimensions  of  the  body.  The  Rayleigh  region  is  characterized  by 
bodies  whose  characteristic  dimensions  are  less  than  a  quarter  of  a 
wavelength.  In  this  region,  the  RCS  is  approximately  proportional  to 
the  volume  squared  and  inversely  proportional  to  the  wavelength  of  the 
fourth  power,  neglecting  degenerate  cases  such  as  discs  and  wires.  The 
resonance  region  covers  the  region  from  approximately  a  quarter 
wavelength  to  a  few  wavelengths  in  body  dimenions.  In  the  resonance 
region  body  shape,  orientation  and  specific  dimensions  are  of 
particular  importance.  The  simplest  example  of  resonance  effects  is 
the  variation  of  cross  section  versus  frequency  curve  for  spheres. 

The  optical  region  covers  bodies  whose  characteristic  dimensions 
are  large  with  respect  to  a  wavelength.  In  the  optical  region  many 
approximation  techniques  such  as  geometric  optics,  physical  optics, 
geometrical  theory  of  diffraction,  and  fringe  wave  theory  are  applied. 

Modeling  of  glint  has  had  less  development  although  the  subject 
has  remained  of  significant  interest  since  the  original  work  by  the 
Naval  Research  Laboratories.  The  subject  of  doppler  scintillation  has 
received  little  attention  although  interest  in  this  phenomenon  has 
increased  in  recent  years  because  of  the  development  of  pulse  doppler 
radar  systems  and  the  increased  resolution  capabilities  desired. 

1.3  Statement  of  the  Problem 

The  evaluation  of  the  performance  of  a  radar  system  requires  a 
knowledge  not  only  of  the  radar  system  hardware  and  characteristics, 


13 


but  also  of  the  radar  target.  Radar  system  evaluations  are  generally 
performed  using  probabilistic  or  statistical  approaches,  with  perfor¬ 
mance  specified  in  terms  like  probability  of  detection,  probability  of 
false  alarm,  and  rms  tracking  accuracy  (in  range,  azimuth  and  elevation 
angles,  and  doppler  frequency).  Such  an  evaluation  requires  statisti¬ 
cal  models  of  the  EM  scattering  parameters  which  affect  system  perfor¬ 
mance.  The  statistical  characteristics  of  the  EM  scattering  parameters 
are,  as  mentioned  in  Section  1.1,  dependent  on  the  target  and  its 
angular  and  translational  motion.  The  objective  of  this  work  is  to 
develop  techniques  to  determine  these  statistics.  In  particular,  this 
work  is  directed  toward  development  of  techniques  to:  (1)  determine 
the  moments  of  the  RCS  and  the  azimuth  and  elevation  errors,  AZER  and 
ELER,  respectively,  and  (2)  determine  the  covariance  functions  for 
these  three  parameters. 

The  work  is  divided  into  two  basic  phases.  The  first  phase, 
described  in  Chapter  2,  is  the  development  of  a  deterministic  model  of 
the  EM  scattering  from  a  target.  The  equations  for  the  RCS  and  target 
induced  error  signals  are  developed  for  an  M-body  target  and  then 
applied  to  an  ellipsoidal  model  of  a  target  drone.  The  second  phase, 
described  in  Chapter  3,  is  the  selection  of  important  statistical 
characteristics  of  the  EM  scattering  and  development  of  techniques  to 
compute  these  statistical  characteristics.  The  characteristic 
function  method  of  computing  these  statistical  characteristics  is 
derived  for  the  discrete  case  and  implemented  on  a  digital  computer. 

A  simulation  was  developed  to  provide  a  check  on  the  general 
accuracy  of  the  analytic  solutions.  A  comparison  of  the  simulations 
and  analytic  solutions  is  presented  in  Chapter  4. 
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A  secondary  objective  of  this  research  is  to  identify  the  areas 


i 


vhere  data  or  theory  is  inadequate  to  support  the  ultimate  objective 
of  developing  stochastic  models  of  the  radar  scattering  for  use  in 
simulations  of  air  defense  systems.  Many  of  these  limitations  are 
Identified  as  they  are  encountered,  and  all  identifiable  problem  areas 
are  covered  in  the  conclusions. 


/pt&x&rtrrn* 
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2.  ANALYTIC  MODELING  OF  RADAR  TARGET  SCATTERING  PARAMETERS 

2.1  General 

In  this  section  the  deterministic  models  for  the  scattering 
parameters  of  interest  are  derived.  These  parameters  are  the  target 
induced  effects  which  cause  errors  in  the  radar  system,  namely  the 
radar  cross  section  (RCS),  the  azimuth  error  (AZER),  the  elevation 
error  (ELER),  and  the  target  phase  (PHAS) .  The  RCS  is  used  here  in 
its  conventional  form.  Glint,  for  the  purposes  of  this  research,  is 
defined  as  the  linear  errors,  in  meters,  in  the  target  azimuth  (AZER) 
and  elevation  (ELER)  positions,  referenced  to  the  target  range.  The 
target  phase  (PHAS)  is  defined  here  as  the  phase  of  the  electromagnetic 
(EM)  vector  reflected  from  the  target  as  compared  to  an  EM  vector 
reflected  from  a  point  scatterer  located  at  the  target  center. 

The  approach  used  assumes  that  the  target  can  be  divided  into 


M  different  scattering  elements  located  at  fixed  points  in  the  target 


coordinate  system.  The  EM  scattering  from  each  element  is  determined 


as  a  function  of  the  aspect  angles.  Then,  the  scattering  parameters 


can  be  determined,  as  functions  of  the  aspect  angles,  by  vector 


•  summation  of  the  scattering  from  the  M  scattering  elements. 

i 

i 

}  Two  primary  coordinate  systems  (shown  in  Figure  2,1)  are  of 

i 

!  interest.  One  is  the  ground  fixed  coordinate  system  with  the  observing 


radar  at  its  origin.  The  azimuth  (A3)  and  elevation  (EL)  angles 


normally  used  by  radar  systems  are  indicated  as  are  the  conventional 
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GROUND  FIXED  COORDINATES 

Figure  2.1.  Coordinate  Systems 
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polar  coordinates.  The  x*  coordinate  is  downrange,  the  y'  coordinate 
is  cross  range,  and  the  z*  coordinate  is  altitude.  The  position  o£  the 
target  is  given  by  (x^,  yQ,  zQ)  or  in  polar  coordinates  by  (Rq,  0q, 

<Pq) .  The  other  coordinate  system  is  the  target  fixed  coordinate  system 
with  the  center  of  the  target  located  at  the  origin.  The  x  axis  is 
the  longitudinal  axis,  the  positive  y  axis  is  out  the  left  wing,  and 
the  positive  z  axis  is  up.  The  yaw,  pitch,  and  roll  angles,  a,  3,  and 
7,  respectively,  are  rotations  about  the  z,  y,  and  x  axes,  respectively 
where  a  and  y  are  defined  in  a  right  sense  and  3  is  in  a  lefthanded 
sense.  The  angles  d  and  <p  which  are  the  polar  coordinates  of  the  radar 
in  the  target  fixed  coordinates  system  are  usually  referred  to  as  the 
aspect  angles. 

The  two  coordinate  systems  are  related  by  the  following 
transformations: 


C11  *12  *13 

*21  *22  t23 

*31  *32  *33 


’x'  -  *o] 

O 

>» 

1 

z*  -  zn 

J 

°J 

*11  *21  t31 
*12  *22  t32 

t13  fc23  *33 


x 

X- 

0 

y 

+ 

y0 

z 

Z- 

J 

m  m 

0 

(2.1) 


(2.2) 


t^  ■  cos  a  cos  3  “  sin  a  sin  3  sin  y 
t^  ■  sin  a  cos  3 

t"13  "  008  a  8*n  ^  +  s*n  a  008  ^  sin  ? 
t2^  ■  -  cos  a  sin  3  sin  y  -  sin  a  cos  3 

C22  "  cos  a  008  ^ 


where 
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t23  «  cos  a  cos  0  sin  7  -  sin  a  sin  0 
t^  a  -  sin  p  cos  7 
t^2  =  -  sin  7 

«  cos  p  cos  7.  (2.3) 

2.2  Scattering  Parameters  of  M-Bodies 

The  modeling  approach  selected  is  based  upon  the  work  of  Dunn  and 
Howard  [28]  and  Ostrovityanov  [24]  .  It  is  assumed  that  there  are  M 
scatterers  on  the  body,  each  of  which  may  be  aspect  angle  dependent. 

The  i-th  scatterer  is  located  at  (x^,  y^,  z^)  in  the  target  coordinate 
system  and  has  RCS  of  S^.  Both  the  location  and  RCS  of  each  scatterer 
are  assumed  to  be  aspect  angle  dependent,  where  the  aspect  angle 
dependence  permits  the  scatterers  to  be  shadowed  or  disappear  from 
view.  Each  scatterer  is  assumed  to  be  independent  of  the  others,  and 
only  single  reflections  are  considered.  This  admittedly  ignores  such 
phenomena  as  interior  corners  and  ducts,  but  these  are  considered  to  be 
of  relatively  minor  importance  statistically.  Polarization  sensitivity 
can  be  included  but  the  derivation  of  the  scattering  equations  assumes 
polarization  independence ,  primarily  for  convenience . 

It  has  been  shown  by  Dunn  and  Howard  [  28J  that  glint  computed  by 
taking  the  normal  to  the  phase  front  is  the  same  as  the  glint  computed 
by  the  ratio  of  the  nonradial  and  the  radial  components  of  reflected 
power.  This  is  true  for  all  currently  used  angle  tracking  systems. 

The  approach  used  here  utilizes  the  Poynting  vector  representation 
of  the  reflected  power.  The  EM  vector  backscattered  from  each  of  the 
M  scatterers  is  computed  and  these  are  summed  vectorially.  The 
Poynting  vector  is  then  calculated  and  decomposed  into  three  orthogonal 
components.  The  first  component  is  radially  directed  from  the  target 
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to  the  radar.  The  other  two  components  are,  for  convenience,  the 
azimuthal  and  elevational  components  in  the  target  coordinate  system. 
The  radial  component  of  the  Poyntlng  vector  is  related  to  the  RCS  and 
the  nonradial  components  are  related  to  the  glint  errors.  The  results 
for  the  scattering  parameters  are  derived  in  Appendix  A  and  are  given 
by: 


M  M 

s  ■  £  £  J  s.  s.  cot  (a.  -  a, \ 
i=>l  >1  1  J  \  1  J/ 

e<p  "  s  £  ^/  Si  S7  f2i  COS  (ai  "  aj) 
€e  “  I  iEi  >i^  Si  SJ  fli  C°S  (ai  ’  aj) 


(2.4) 


(2.5) 


(2.6) 


PHAS  «*  arc  tan 


/  £  J  S.  sin  a, 
!  1=1 


«  ______ 

£  /  s,  cos  a 
j»l  J  J 


(2.7) 


where  and  e  are  the  errors  in  the  target  coordinate  system,  0  and 

y  p 

cp  are  the  target  aspect  angles,  and 

«  •—  ^-x^  sin  0  co3  cp  -  sin  0  sin  cp  -  cos  0^  (2.8) 

fj^  ■  -x^  cos  0  cos  cp  -  y^  cos  0  sin  q>  -  sin  0  (2.9) 

f2£  ■  -x^  sin  cp  +  y^  cos  cp  •  (2.10) 

The  €  and  e  errors  given  in  (2.5)  and  (2.6)  must  be  transformed 

<p  0 

to  azimuth  and  elevation  errors,  AZER  and  ELER  respectively,  in  radar 
coordinates  to  account  for  the  orientation  of  the  target.  Let 
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[D(0,<p)l 


sin  9  cos  cp  sin  9  sin  cp 
cos  9  cos  cp  cos  9  sin  cp 


-  sin  cp 


cos  cp 


cos  9 
•sin  9 
0 


Then,  the  azimuth  and  elevation  errors  are  obtained  £rom 


€ 

0 

r  /  \i 

r  -i  i 

-i  i 

ELER 

“  [D  (V  ,Po)J 

b  1 

d  (e,  cp) 

ee 

AZER 

e 

cp 

(2.11) 


(2.12) 


€  is  a  dummy  variable,  and  [T  ]  is  [T  ]  with  the  mean  values  of  a, 
P,  and  y  substituted  for  a,  p,  and  y .  For  3traight  and  level  flight 
where  the  mean  pitch  and  roll  angles  are  zero 


and  (2.12)  yields 


AZER  =»  € 

<P 


ELER  =s  €0 


(2.13) 


In  order  to  derive  the  and  (x^  y^,  z^)  of  each  scattering 
element,  the  target  is  modeled  by  M  ellipsoids. 

2.3  Ellipsoidal  Surface  Model 

It  is  assumed  that  the  surfaces  of  the  target  can  be  approximated 
by  ellipsoids,  as  depicted  in  Figures  2.2  and  2.3  for  a  target  drone 
BQM-34A.  Ellipsoids  were  chosen  because  of  the  relative  ease  with 
which  one  can  find  and  describe  the  scattering  from  the  specular  point. 
Each  ellipsoid  is  described  by  its  conjugate  radii  and  its  position 
and  orientation  in  the  target  fixed  coordinate  system.  An  associated 
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modulating  function  is  defined  for  each  ellipsoid  which  truncates  the 
ellipsoid,  approximates  the  shadowing  by  the  other  surfaces,  and 
provides  a  relatively  smooth  transition  into  the  field  of  view.  It 
may  also  be  used  to  provide  an  empirical  method  of  improving  the  fit 
of  calculated  to  empirical  data.  Each  ellipsoid  may  be  considered  as 
a  point  scatterer  with  RCS,  S^,  and  location  (x^,  y^,  z^)t  determined 
by  the  specular  point  [41,  42] . 

First  a  single  ellipsoid  located  at  the  origin  with  its  principal 
axes  oriented  along  the  Cartesian  coordinate  axes  is  considered.  The 
RCS  of  such  an  ellipsoid  is  given  approximately  by 


S  «  it  R1  R2 


(2.14) 


where  Rj^  and  R£  are  the  principal  radii  of  curvature  at  the  specular 
point  [41].  The  specular  point  is  found  by  finding  the  point  on  the 
surface  of  the  ellipsoid  having  the  same  direction  cosines  as  the 
lines  of  sight. 

Let 


2  2  2 

F(x,  y,  z)=^  +  ^  +  ^-  l  =  0  (2.15) 

a  b4  c 


be  the  equation  of  an  ellipsoid.  Then  the  normal  to  the  ellipsoid  at 
(x,  y,  z)  is  given  by 


so  that  the  direction  cosines  of  the  normal  are  given  by 


cos  5 
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X 

2 

a  r 


(2.16) 
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(2.17) 


However,  the  direction  cosines  for  the  line  of  sight  are  given  by 


cos  5  =  sin  9  cos  cd 
x  Y 

cos  &  **  sin  9  sin  m 

y  v 

cos  &  ■  cos  6  (2.18) 

z 

where  9  and  cp  are  the  aspect  angles  of  the  ellipsoid  given  in  its  own 

(local)  coordinate  system.  The  ambiguity  of  the  normal  can  be  resolved 

by  noting  that  the  outward  normal  must  be  in  the  same  octant  as  the 

line-of-sight.  Thus  from  (2.17)  and  (2.18)  the  parametric  angular 

coordinates,  u  and  v  ,  of  the  specular  point  are  obtained  as 
s  s 

v  *  arctan !  7  tan  cp 

8  l  cL 

1  f~2  2  2  2  l 

u  =>  arctan  —  tan  9  J  a  cos  cp  +  b  sin  cp  /  .  (2.19) 

s  |c  j 

The  product  of  the  principal  radii  of  curvature,  X  K, ,  is 

derived  as  in  [41]  by 

specular 

point  (2.20) 

where 
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The  application  of  (2.20)  to  the  ellipsoid  (2.15)  results  in  the 
expression 


_  _  2.2  2(X3±ys1zs 

h  h  "  a  b  0  +  +  - 

\a  b  c 


(2.22) 


where  (xg,  yg,  zg)  are  the  coordinates  of  the  specular  point  and  are 
given  from  (2.16)  by: 


x  =  a  sin  u  cos  v  ■  —  sin  0  cos  m 
s  s  s  p  v 


y  ■  b  sin  u  sin  v  *  —  sin  0  sin  rp 
•'s  s  3  P  v 


z  =  c  cos  u  *  —  cos  0 
s  s  p 


(2.23) 


where 

y  22  2  22  2  22 

a  sin  0  cos  cp  +  b  sin  0  sin  cp  +  c  cos  0 

Consequently,  the  RCS  of  the  ellipsoid  may  be  found  from  (2.14)  and 
(2.22)  to  be  given  explicitly  by 

„  2  .2  2 

S  =  *  Rx  R2  *  ■  ■  a'  |  C  .  (2.24) 

P 
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The  results  given  in  (2.23)  and  (2.24)  for  a  single  ellipsoid  may  now 
be  applied  to  each  of  the  M  ellipsoids,  after  the  proper  coordinate 
trans  formations . 

Consider  now  an  ellipsoid  which  has  its  center  located  at 
(*ti»  yti,  zt^)  in  the  target  coordinate  system.  A  local  coordinate 
system  for  the  ellipsoid  is  defined,  with  the  origin  located  at  the 
center  of  the  ellipsoid  and  with  the  coordinate  axes  coincident  with 
the  principal  axes  of  the  ellipsoid.  .  Under  these  conditions  the 
ellipsoid  is  given  by  (2.15).  The  angular  rotations  from  the  target 
coordinate  system  to  the  local  coordinate  system  are  defined  as  0fc  and 
(pt,  where  cpfc  is  taken  first  about  the  z  axis  and  then  0 is  taken  about 
the  new  y  axis. 

Let  x,  y,  and  z  be  the  coordinates  of  a  point  in  the  target 
fixed  coordinate  system,  and  let  x' ,  y' ,  and  z'  be  the  coordinates  of 
the  point  in  the  local  coordinate  system.  Then  the  coordinates  of  the 
point  are  related  by  the  following  transformation: 


1  _ 

1 _ 

X  “  Xti 

y' 

«  [TJ 

y  ■  yti 

z’ 

m  m 

z  -  zti_ 

where 


[TJ 


cos  <j>t  cos  et  sin  cpt  cos  -sin  0t 

-sin  cpt  cos  cp  0 

cos  (pt  sin  0t  sin  cpt  sin  cos  0fc 


(2.25) 


Let  0  and  cp  be  the  aspect  angles  of  the  target  and  (x^,  y^,  z^) 
be  the  coordinates  of  the  point  on  the  unit  sphere  with  these  angular 
coordinates,  then 
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If  r,  s,  and  t  are  the  coordinates  of  this  point  in  the  local  coordi- 
nate  system  of  the  ellipsoid,  then 
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sin  0  sin  cp 
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cos  0 
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(2.27) 


hence  the  local  aspect  angles  can  be  found  from 


cpt  =  arctan^-^ 

0^  »  arctan^-'^/r2  +  s2  ^  .  (2.28) 

The  octant  in  which  the  angles  are  located  must  be  determined  by 
consideration  of  the  signs  of  r,  s,  and  t. 

2.4  Computational  Procedure 

Let  0  and  cp  be  the  target  aspect  angles;  0 ^  and  cp^  be  the  local 
aspect  angles  of  the  i-th  ellipsoid;  0fci  and  cpt^  be  the  angular  trans¬ 
lations  from  the  target  coordinate  system  to  the  local  coordinate 
system  for  the  i-th  ellipsoid;  (x^,  yt^»  z^)  be  the  location  of  the 
center  of  the  i-th  ellipsoid;  (x|,  y^,  z'^)  be  the  coordinates  of  the 
specular  point  in  the  local  coordinate  system  of  the  i-th  ellipsoid; 
and  x^,  y^,  and  z^  be  the  coordinates  of  the  specular  point  of  the 
i-th  ellipsoid  in  the  target  fixed  coordinate  system. 

The  starting  point  for  each  computation  is  the  set  of  target 
aspect  angles  0  and  cp.  For  each  of  the  N  ellipsoids,  the  following 
procedure  is  used: 


(1)  Transform  from  target  aspect  angles  to  local  aspect  angles 
using  (2.28). 

(2)  Determine  the  location  of  the  specular  point  in  the  local 
coordinate  system  using  (2.23). 

(3)  Determine  the  RCS,  Sj,  of  the  specular  point  using  (2.24). 

(4)  Transform  the  position  of  the  specular  point  from  local 
coordinates  (x^,  y^,  z|)  to  target  coordinates  (x^  y^  z^)  by 


(5)  Determine  the  magnitude  of  the  modulating  function  G(i)  and 
the  value  of  the  RCS  to  be  used  in  (2.4)  through  (2.7)  by 


St  *  G(i) 

The  values  of  the  scattering  parameters  for  the  aspect  angles  6 
and  <p  are  found  by  applying  (2.4)  through  (2.7)  where  the  location 
(*i,  y^,  zt),  and  RCS,  S^,  of  the  i-th  scatterer  are  determined  as 
described  above. 

2 , 5  Numer ica 1  Examp 1 e 

The  results  of  the  last  three  sections  will  now  be  applied  to 
the  special  case  of  the  BQM-34A  Target  Drone.  Figures  2.2  and  2.3 
depict  the  actual  Target  Drone  and  the  ellipsoidal  approximation.  The 
ellipsoidal  model  in  this  case  consists  of  eight  ellipsoids.  The 
configuration  (as  defined  by  the  conjugate  radii),  location  (xt,  yfc, 
zt)  and  orientation  (6t,  cpt)  for  each  ellipsoid  is  given  in  Table  2.1. 
The  modulating  functions  for  these  ellipsoids  are  given  in  Table  2.2. 


so 


The  procedures  derived  in  Sections  2.2  and  2.4  are  used  to  compute 
the  scattering  parameters  as  functions  of  tb~  aspect  angles,  6  and  cp. 
Typical  results  are  shown  in  Figures  2.4,  2.5,  and  2.6. 

Figure  2.4  is  the  calculated  RCS  for  L-Band  plotted  in  decibels 
referenced  to  l-square  meter  (RCSdBM)  versus  yaw  angle  for  zero  roll 
and  zero  pitch.  Zero  yaw  angle  corresponds  to  nose-on,  90  degrees 
corresponds  to  broadside  and  180  degrees  corresponds  to  tail-on.  As 
compared  to  measured  data,  the  calculated  RCS  in  the  angular  regions 
from  0  to  about  35  degrees  and  from  about  140  to  180  degrees  is  low. 

In  the  headon  region  the  discrepancy  is  about  15  decibels  and  in  the 
tail-on  region  about  10  decibels.  These  discrepancies  are  due  largely 
to  two  factors.  The  first  factor  is  the  simplicity  of  the  modeling 
approach  used.  The  second  factor  is  the  transparency  of  certain  parts 
of  the  target  drone.  The  analytic  model  ignores  the  engine  intake  and 
exhaust  ducts  thereby  ignoring  significant  scattering  over  certain 
aspect  angles.  The  data  on  which  the  comparison  is  based  were  measured 
on  a  target  drone  with  considerable  transparent  portions.  In  addition 
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to  the  radome,  the  front  cowl  forward  of  the  engine  is  transparent. 
Therefore,  the  RCS  within  about  35  degrees  of  nose-on  is  the  RCS  of  the 
bulkheads,  forward  engine  structure,  autopilot  and  electronics  packages, 
and  other  components.  The  end  of  the  tail  cone  which  houses  the 
recovery  parachute  is  also  transparent.  These  factors  make  a  comparison 
extremely  difficult.  However,  the  lobing  structure  is  approximately 
correct  as  is  the  broadside  RCS. 

Figures  2.5  and  2.6  are  the  azimuth  and  elevation  errors,  AZER 
and  ELER  respectively,  expressed  as  the  linear  error  at  the  target. 

The  conditions  are  identical  to  those  for  Figure  2.4. 


<p  (deg) 


Figure  2.5.  Azimuth  Error  e  Versus  Yaw  Angle  cp 


€q  (meters) 


V  (deg) 


Figure  2.6. 


Elevation  Error  e 

e 


Versus  Yaw  Angle  cp 
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3.  THE  STATISTICAL  PROPERTIES  OF  THE  SCATTERING  PARAMETERS 

3.1  General 

The  primary  objective  of  this  research  has  been  to  develop  tech¬ 
niques  for  computing  the  statistical  parameters  of  a  radar  target,  in 
particular  an  aircraft  type  target.  These  data  can  be  used  to 
generate  a  vector  stochastic  process  representing  the  radar  scattering 
parameters  of  the  target  which  can  be  used  in  system  simulations  and 
evaluations.  The  data  may  also  be  used  in  the  design  of  optimum 
filters  for  tracking  systems. 

Section  3.2  derives  the  statistics  of  the  aspect  angles  using 
approximation  techniques.  Section  3.3  applies  the  aspect  angle 
statistics  to  the  analytic  scattering  model  developed  in  Chapter  2 
to  derive  the  statistical  properties  of  the  scattering  parameters. 

3.1.1  The  Statistical  Properties  Under  Investigation 

The  available  techniques  for  obtaining  a  dynamical  repre¬ 
sentation  of  a  stochastic  process  require  that  the  probability  density 
function  and  autocorrelation  function  be  known,  and  that  these  functions 
satisfy  certain  criteria.  These  techniques  can  be  generalized  to  N 
dimensions,  when  the  matrix  of  covariance  functions  is  known.  One 
general  approach  requires  only  that  the  probability  density  function 
have  finite  variance  and  that  the  covariance  functions  be  written  as 


sums  of  exponentials. 


The  RCS,  AZER,  and  ELER  are  complicated  nonlinear  functions  of 
the  aspect  angles  as  indicated  in  Chapter  2  so  that  the  conventional 
transformation  of  variables  approach  to  obtain  the  probability  density 
function  is  impractical  if  not  impossible.  Therefore)  other  properties 
of  the  probability  density  functions  must  be  investigated.  The  most 
descriptive  properties  which  can  be  determined  are  the  moments.  A 
complete  description  of  a  probability  density  function,  in  terms  of  its 
moments,  requires  that  all  of  the  moments  be  known.  For  a  Gauss iar 
density,  all  of  the  moments  can  be  expressed  in  terms  of  the  mean  and 
variance.  These  two  moments  are  not  sufficient,  in  general,  to 
characterize  an  unknown  density  function.  However,  computation  of  all 
moments  is  not  a  practical  undertaking  so  some  compromise  must  be  made. 
A  reasonable  compromise  appears  to  be  the  first  four  moments.  The 
third  moment  measures  the  skewness  or  asymmetry  and  the  fourth  moment 
measures  the  shape  of  the  distribution  as  compared  to  the  Gaussian 
distribution. 

These  moments  can  be  used  to  approximate  the  fit  of  an  unknown 
probability  density  function  by  use  of  known  techniques  including  the 
use  of  Gram-Charlier  series  [57]  where  the  unknown  density  functions 
do  not  differ  markedly  from  Gaussian  density  functions  or  those 
derivable  from  Gaussian  density  functions.  Another  technique  is 
described  by  Wragg  and  Dowson  [58]  for  density  functions  on  [0,  «>), 

The  statistical  properties  of  the  scattering  parameters  under 
investigation,  then  are  the  first  four  moments  and  the  matrix  of 
covariance  functions.  Since  it  is  not  possible  to  assign  probabilis¬ 
tic  properties  to  a  target  trajectory,  the  statistics  to  be  derived  are 
all  for  a  given  trajectory.  Therefore,  the  moments  and  the  correlation 
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functions  are  defined  as  conditional  expectations  given  the  flight 
path  in  the  ground  fixed  coordinate  system. 

3.1.2  The  General  Approach 

The  time  history  of  the  scattering  parameters  of  an  air¬ 
craft  cannot  be  represented  by  a  stationary  process.  The  means  and 
variances  of  the  scattering  parameters  are  time-varying  functions  of 
the  aspect  angles  0  and  <p,  which  are  in  general  functions  of  the 
trajectory  and  hence  of  time.  The  relationship  among  the  aspect  angles 
and  the  target  position  and  attitude  in  the  ground  fixed  coordinate 
system  .ire  given  by 

••  «« 
sin  0(t)  cos  <p(t) 

sin  0(t)  sin  cp(t) 

cos  0(t) 

where  [T]  is  given  in  (2.3). 

The  target  coordinate  angles  0g(t)  and  rp^(t)  describe  the  target 
position  in  space  as  measured  in  the  ground  fixed  coordinate  system. 

The  Euler  angles,  a,  P,  and  7,  appearing  in  [T]  are  functions  of  the 
target  motion,  making  the  [T]  a  function  of  time  for  a  known 
trajectory.  Let  a(t),  f3(t),  and  t)  be  the  Euler  angles  corresponding 
to  a  given  target  trajectory  (position,  velocity,  and  acceleration). 

If  these  are  substituted  in  [T]  in  (3.1)  the  resulting  0(t)  and  (p(t) 
are  then  purely  deterministic  functions  of  time.  However,  in  general, 
the  Euler  angles  a,  p,  and  7  contain  a  stochastic  component  due  to 
perturbations  in  the  motion,  wind,  vibrations,  etc.  For  simplicity 
it  may  be  assumed  that..-the  stochastic  components  of  a,  p,  and  7  are 
additive  with  zero  mean  and  exponential  autocorrelation  function. 


-sin  0Q(t)  cos  cpQ(t) 

=  [T(t)J 

-sin  0 q ( t )  sin  cpQ(t) 

-cos  0Q(t) 

(3.1) 


i 

l 


'«*  J*?  «  » i*-»J*.  ***',  .  *»-.  - 
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Therefore,  a(t),  (3(t),  and  y(t)  are  assumed  to  have  time  varying 
means  3(t),  £(t),  7(t)  and  autocorrelation  functions. 


Rv(tl’ 


c*)  ■  7('i)  ~(S) 


+  a 


(3.2) 


v  -  a,  p,  7 


The  perturbations  have  generally  small  time  constants  as  compared  to 
the  time  constants  of  aircraft  trajectory.  Consequently,  a(t),  p(t), 
and  7(t)  may  be  considered  to  be  wide-sense  stationary  for  |tg  -  t- 1 
small  enough.  Now  d(t)  and  <p(t)  in  (3.1)  are  given  by  time  varying 
nonlinear  transformations  of  a,  p,  and  y.  Since  the  time  constants  of 
the  stochastic  components  are  small  compared  to  the  time  constants  of 
the  trajectory,  and  hence  of  the  transformation,  it  can  be  shown  that 
0(t)  and  cp(t)  are  therefore  time  varying  random  processes  whose  means 
and  autocorrelation  functions  have  the  form 


E(0(t)}  =  SXt)  *  0(t) 

E{<p(t)}  ■  "cp(t)  35  q5(t)  (3.3) 

R*  (‘i.  *2)  *  W)  +  ",  h) 

R(p  (ei*  h)’  v(tlfS(^2)  *  0l(tl)°ri;(t2),;5(t2  ‘  Cl)  • 

Here,  it  has  also  been  a.,  ^ed  that  the  stochastic  perturbation  of  the 
Euler  angles  are  relatively  small.  Consequently  for  a  given  flight 
path,  Q (t)  and  <p(t)  have  time  varying  means  and  variances.  Further¬ 
more,  if  } ^2  -  tj |  is  small,  then  their  covariance  is  approximately 
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'«('!•  t2)“  ‘’e(ti)»(t2  '  ci) 

h)~  fyifo  ~  *i)  <3,5) 

while  for  large  |t2  -  they  are  uncorrelated.  Similarly,  for 
|t2  -  tjJ  small  enough  that  p(t2  -  t^)  is  not  negligible,  0(t^  =  0(t2) 
and  *  cp(t2).  Hence,  they  are  wide-sense  stationary  for  suffi¬ 

ciently  small  time  intervals. 

The  approach  to  estimating  the  aspect  angle  statistics  is  to 
expand  the  implicit  expressions  for  0  and  cp  contained  in  (3.1)  about 
the  deterministic  values  0  and  cp  which  are  obtained  by  substituting 
a,  3,  and  y  into  (3.1)  for  a,  3,  and  y.  The  order  of  the  approxima¬ 
tion  will  be  reduced  to  the  minimum  acceptable  level. 

3.2  Aspect  Angle  Statistics 

The  derivation  of  the  statistical  characteristics  of  the  aspect 
angles,  Q  and  cr»  requires  a  knowledge  of  the  target  position  in  the 
ground  fixed  coordinate  system  and  the  target  motion  in  each  of  its 
6  degrees  of  freedom.  The  translational  motion  is  specified  in  the 
ground  fixed  coordinate  system.  The  angular  motion  is  more  difficult 
to  define  but  relates  the  angular  position  of  the  target  fixed 
coordinate  system  to  the  ground  fixed  coordinate  system.  Neglecting 
any  random  angular  motions  it  is  assumed  that  the  target  flies  with 
its  longitudinal  axis,  x  axis,  coincident  with  the  velocity  vector 
and  its  vertical  axis,  z  axis,  opposing  the  combined  acceleration  due 
to  gravity  and  any  maneuvers  other  than  longitudinal  acceleration. 
(This  approach  clearly  ignores  aerodynamic  theories  related  to  flight 
[59]  but  permits  a  reasonable  analysis  without  requiring  a  detailed 
aerodynamic  model  of  the  target.)  The  yaw,  pitch,  and  roll  angles, 
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a,  p,  and  y  respectively,  required  to  go  from  axes  parallel  to  the 
ground  fixed  Cartesian  coordinate  axes  to  the  target  fixed  Cartesian 
coordinate  axes  are  shown  in  Figure  3.1. 

Let  x',  y',  and  z'  be  the  Cartesian  coordinate  axes  of  the  ground 
fixed  coordinate  system  and  x,  y,  and  z  be  the  Cartesian  coordinates 
of  the  target  fixed  coordinate  system.  As  defined  previously,  (1)  the 
x  axis  is  the  longitudinal  or  roll  axis  of  the  target  with  the  velocity 
vector  pointing  in  the  positive  direction;  (2)  the  y  axis  is  the  pitch 
axis;  and  (3)  the  z  axis  is  the  yaw  axis.  The  yaw,  pitch,  and  roll 
angles  (Euler  angles)  ,  a,  p,  and  y  respectively,  are  defined  such  that 
a  and  y  are  conventional  right-handed  angles,  but  p  is  left-handed. 

Let  x*,  y*,  and  z*  be  the  intermediate  axes  during  the  sequence  of 
angular  transformations  from  one  coordinate  system  to  the  other.  The 
transformation  between  coordinate  systems  is  unique  only  when  the  order 
of  the  angular  transformations  is  defined.  In  this  work,  the  order 
is  defined  to  be  p,  y,  and  a,  in  that  order.  Figure  3.1  depicts  the 
angular  coordinate  transformation  from  the  ground  fixed  coordinates  to 
the  target  fixed  coordinates.  The  first  rotation  is  p  about  the  y' 
axis,  the  second  rotation  is  y  about  x*,  the  intermediate  x  axis,  and 
the  third  is  a  about  the  final  z  axis. 

The  aspect  angles  are  related  to  the  target  position  and  motion  by 
(3.1)  where  [T]  is  defined  in  (2.3).  If  the  a  rotation  is  factored 
out  of  <p,  it  is  noted  that  a  is  additive  to  cp  resulting  in  a  simplifi- 
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It  is  possible  now  to  find  an  exact  analytic  solution  for  the 
joint  probability  density  function  for  9  and  cp»  as  given  in  Appendix 
B.  Clearly  the  equations  in  the  right  hand  matrix  are  nonlinear  and 
do  not  permit  a  visualization  of  the  results  even  when  the  density 
functions  for  a,  p,  and  y  are  given.  Numerical  results  can  be  obtained 
by  implementing  standard  transformation  of  variables  techniques  on  a 
digital  computer.  An  even  more  formidable  limitation  exists.  Although 
it  is  possible  to  obtain  the  moments  by  numerical  integration,  it  is 
not  possible  to  obtain  the  covariance  functions  for  9  and  cp  which  are 
necessary  in  computing  the  covariance  functions  of  the  scattering 
parameters . 

The  aspect  angles  9  and  cp  can  be  expanded  in  Taylor  series  about 
their  quas i -means  #(t)  and  'cp(t).  For  notational  convenience,  the 
time  dependence  of  the  means  and  variances  is  being  dropped  from  the 
notation  for  the  rest  of  the  analysis.  If  variances  of  the  random 
components  of  the  Euler  angles  are  sufficiently  small,  9  and  cp  can 
be  approximated  by  the  first  two  terms  of  the  Taylor  series  expansions. 


(3.7) 
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be 

.  be 

A 

dp 

7s  ^7 

+  Vs  W 

P.7 

P.7 

P.7 


i  „2  re 
dp2 


+  2  *s  -n„2 


P.7 


1  2  o0 


P.7 


(3.8) 


where 


a  =  a  -  a 
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P3  "  P  -  P 


7.  =  7  “7 


By  use  of  implicit  differentiation  and  algebraic  manipulations  the 
partial  derivatives  can  be  shown  to  be: 


be 

5p 


__  _  =  COS  7  cos  Cp 

P.7 


(3.9) 


d2e 

dp2 


1  ~  2.  “  ~  2  ~ 
*  j  sin  2/  sin  cp  +  j  (1  +  cos  2/)  cot  0  sin  cp 


P.7 


(3.10) 


be 

57 


sin  cp 


(3.11) 


P.7 


2 

d  0  I  ~  2  ~ 

=  cot  0  cos  cp 


^7 


(3.12) 


P.7 


a2e 

5p57 


'V  »v 


-  cos  cp  (sin  7  -  cos  7  cot  0  sin  cp)  (3.13) 


P.7 
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rv 


-  sin  7  -  cos  7  cot  0  sin  cp  (3.14) 


P.7 


e.r 


2  sin  ly  cot  9  cos  9  +  T  sin  2tp(l  +  cos  2y)(l  +  2  cot  0) 


=  cos  cp  cot  0 


(3.15) 

(3.16) 


P.7 


£J2 

^72 


=  -  j  sin  2cp  (1  +  2  cot2  0) 


(3.17) 


P.7 


f\*  <v 


9  ^  0  ^  rvi 

*  sin  7  sin  cp  cot  0  -  cos  y  (cos  <p  +  cot  0  cos  2  cp) 

(3.18) 

The  accuracy  of  these  expansions  may  be  examined  by  obtaining  the 
statistics  of  the  resulting  series: 


0»E(0]«0+|op  £sin  2  7  sin  %  +  ^1  +  cos  2  yj  cot  0  sin2 


2  ~  2  ^ 
+  cot  0  cos  cp 


5] 

(3.19) 


1  2  r 

2  L 


V  =  E[cp}  55  cp  +  7  Oq  I  sin  2  y  cot  0  cos  cp 


j  sin  2  cp  ^1  +  cos  2  7  j  ^1  +  2  cot  0^| 


12  «  ~  /  2  ~\ 
j  sin  2  cp  ^1  +  2  cot  0J 


(3.20) 
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o 2  m  e|(0-?)2|  =  0^  cos2  y  cos2  ©  +  o2  sin2  cp 


e 


22  2  ~  (  —  —  ~ 

+  0^  o^  cos  cp  ^sin  y  -  cos  y  cot  0  sin 


?y 


5] 


1  4  r  —  /v  ^  'w  2  ~ 

jr  Op  j^sin  2  7  sin  ©  +  (1  +  cos  27)  cot  0  sin 

4  2  ~  2  ~ 

+  0^  cot  0  COS  cp 


2  (/  — \  21  2  2  —  —  ~  ~  2 

°cp  53  E\v  “  fa  aa  +  ap  (8*n  ^  +  cos  ?  cot  0  s*n 


(3.21) 


2  2  ~  2~  2  2  I" 

+  0^  cot  0  cos  cp  +  Op  0^  |^sin  7  sin  cp  cot  0 
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—  «  2  ~  2  ~  \1 
cos  7  (cos  cp  +  cot  0  cos  2  cp  I 


+  £  Op  Jsin  2  7  cot  0  cos  cp  +  j  sin  2  ©  (1  +  cos  2  7)^1  +  2  cot2  0 


+  ^  sin2  2  ©  ^1  +  2  cot2  0^ 


r 

(3.22) 


30°cp  p0,cr  a  E^0  "  0)(«p  “  Cp)}  =  Op  cos  7  cos  cp 
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*N/  ~ 


X  (-  sin  7  -  cos  7  cot  0  sin  co)  +  j  0  sin  2  ©  cot  0 


I  L[  —  «—  ~ 

+  £  Op  j^sin  2  7  sin  ©  +  (1  +  cos  27)  cot  0  sin  cp  J 

|sin2  7  cot  0  cos  ©  +  |,sin2  ©  (l  +  cos  2  7)(1  +  2  cot2  0)j 

2  2  «*~  ••  ^ 

-  OpO^  cos  cp  (sin  7  -  cos  7  cot  0  sin  cp) 

[—  ~  —  /  2  ~  2  ~  ^  \*1 

sin  7  sin  ©  cot  0  -  cos  7  lcos  cp  +  cot  0  cos  2  cp  I 
lA  ~  2  **  ~  /  2  ^  \ 

•  j  a  cot  0  cos  cp  sin  2  cp  ^1  +  2  cot  6  j 


(3.23) 
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Some  knowledge  of  the  statistics  of  the  Euler  angles  is  necessary 
in  order  to  establish  bounds  on  the  errors,  and  hence  evaluate  the 
accuracy  of  the  expansions. 

There  is,  unfortunately,  a  lack  of  published  or  available  data 
on  the  angular  motions  of  an  aircraft  in  flight.  An  aircraft  is 
designed  to  fly  and  perform  a  mission  safely.  Only  those  flight 
characteristics  which  affect  safety  or  the  ability  to  perform  the 
specified  mission  are  investigated  and  reported.  The  small  scale 
random  motions  of  interest  here  are  generally  ignored.  The  pilot,  or 
autopilot,  can  handle  these  as  a  minor  part  of  the  normal  flight 
activities.  The  coupling  among  yaw,  pitch  and  roll,  especially  when 
flying  with  the  autopilot  operating,  is  familiar  to  many.  It  would 
seem  reasonable,  based  upon  this  and  other  factors,  to  expect  that  the 
quasi-random  angular  motions  would  be  similarly  coupled.  However,  the 
lack  of  data  make  it  impossible  to  determine  what  the  coupling  or 
covariance  function  should  be. 

It  is  necessary  to  assume  some  probability  and  covariance 
functions  for  the  Euler  angles.  The  most  useful  family  of  probability 
density  functions  for  angles  is  that  used  by  Viterbi  [60] .  It  is 
given  by 

-  eX2-(X0(s)  B)  •  0S9S2n  (3.24) 

where  6  is  the  angle  and  s  is  a  parameter  related  to  the  variance  of 
the  distribution.  For  s  =  0,  the  density  is  uniform  over  the  range 
of  -it  to  t,  and  for  s  =>  »,  the  density  is  a  delta  function  at  Q  *  0. 
This  density  is  closely  approximated  by  a  normal  density  for  s  >>  1, 
as  shown  by  Viterbi.  Due  to  the  lack  of  published  data,  it  is  assumed 


m 
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that  the  random  angular  motions  of  the  target  are  independent  normal 
(Gaussian)  processes  with  exponential  autocorrelation  functions.  This 
is  an  arbitrary  and  perhaps  simplistic  approach,  but  represents  a 
practical  and  mathematically  tractable  approach.  The  choice  of  a 
normal  process  is  based  upon  the  fact  that  the  random  motions  are 
small.  The  only  justification  for  the  exponential  autocorrelation 
function  is  mathematical  convenience. 

For  a  basically  stable  aircraft  under  cruising  flight  conditions, 
standard  deviations  of  the  Euler  angles  of  less  than  100  milliradians 
or  about  5.7  degrees  seem  reasonable.  Autocorrelation  time  constants 
on  the  order  of  1.0  to  5.0  seconds  also  seem  reasonable. 

Without  loss  of  generality,  but  in  the  interest  of  simplifying  the 
algebra,  it  is  assumed  that  the  aircraft  is  flying  a  rectilinear  flight 
path.  This  implies  that  p  =  7  =  0,  and  that  0  =  0Q  and  9  =>  cpQ  +  a. 

From  practical  considerations,  the  area  of  prime  interest  in  the 
elevational  aspect  angle  0  is  for  n/2  £  0  <  3it/4.  In  this  region 
-1  £  cot  0  £  0.  With  this  information  it  is  possible  to  evaluate  the 
general  accuracy  of  the  Taylor  series  expansions  of  0  and  c',  (3.7) 
and  (3.8).  Equations  (3.19)  through  (3.23)  can  be  written  as: 


0  =  0  + 


1  ~  /  2  2~  2  2~\ 

2  co  .  0  I  Og  sin  cp  +  0^  cos  cp  j 

cp  =  cp  +  ^  sin  2  rp  ^1  +  2  cot2  0^j2  "  °7  ^ 

2  2  2  ~  2  2~122  2~  2~ 
°0  "  °P  C°S  ^  +  °7  3in  9  +  4  COt  0  Sitl  ^  ? 


2  ~  /  4  4^  4  4  ~ 

+  cot  0  I  Op  sin  cp  +  0^  cos  cp 


) 


(3.19A) 

(3.20A) 


(3.21A) 


47 


2  2,  2-/2  .  2  ~  2  2  — \ 
j  ■  +  cot  0  sin  cp  +  Oy  cos  cp  J 


22/  2  ~  2- 
+  ^cos  cp  +  cot  0  cos  2 


’) 


2 

+  «|  sin2  2  qj  (l  +  2  cot2  e)  (o4  +  a4)  (3.22A) 


•%*  <v 


0.  o  p*  »  ■»  sin  2  cp  cot  0 

e  ye, cp  2  ^ 


PI  V 

7 ' 

(°7  °p) 


q)  ^1  +  2  cot2  0^( 


•w  -  —  ^  a  2  4  2 

+  -f  cot  0  sin  2  cp  (1  +  2  cot“  0  j(o^  sin  <p  -  o^  cos 


?) 


-  ^  o2o2  sin  2®  cot  0  ^cos2  cp  +  cot2  0  cos  2  cp^  .  (3.23A) 


By  inspection,  it  is  clear  that  the  second  order  terms  generally  add 

less  than  0.1  sup  ( o ,  a  )  to  the  means  and  standard  deviations  of  the 

P  7 

aspect  angles.  First  order  (linear)  approximations  will  be  sufficient 
for  this  analysis.  The  following  equations  will  be  used  in  the 
computations  later  in  this  chapter. 


0  *  0n  +  B  cos  q)  +  7  sin  co  (3.25) 

US  s 

cu  *  *'q  +  a  -  cot  0q  sin  co  +  7g  cot  0q  cos  q)  (3.26) 

where  co  ■  co^  +  a  ,  The  means  and  covariance  functions  of  the  aspect 
angles  are  given  by: 


0  -  0Q 

(3.27) 

cn  =  cOq  +  a 

(3.28) 

S^l’  fc2 

)  ‘°l 

2  "kBT  2  2  "kyT 

cos  cp  e  p  +  sin  q5  e  ' 

(3.29) 

\  2 
,  i  “0 

/  °a 

-k  t  2  2  2  -V 

e  +  0^  cot  0q  sin  q>  e  p 

,2  2  2  _  "  V 

+  cot  0q  cos  cp  e 


(3.30) 
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C8„  (*!•  *2)  -  lsi“2  *  e0  \4  -  4  e'kpT)  (3>31) 

where  t  =  I  fc2  "  t]_l  and  ^a»  and  are  the  inverses  of  the  time 
constants  of  the  autocorrelation  functions  of  a,  f},  and  7, 
respectively.  The  linearization  also  implies  that  the  joint  density 
function  of  e  and  <p  is  approximately  Gaussian. 

3.3  Statistics  of  the  Scattering  Parameters 

3.3.1  Derivation  of  General  Approach 

The  previous  chapter  and  previous  sections  of  this  chapter 
have  developed  the  necessary  background  and  inputs  to  the  derivation  of 
the  statistics  of  the  scattering  parameters.  The  following  development 
holds  true  for  all  of  the  scattering  parameters,  so  only  the  RCS  will 
be  mentioned.  The  application  to  the  other  scattering  parameters  will 
be  made  when  appropriate. 

As  mentioned  in  Section  3.1,  it  is  not  possible  to  obtain  the 
probability  density  function  for  the  RCS  using  conventional  transfor¬ 
mation  of  variables  techniques.  Therefore,  the  first  four  moments  of 
the  probability  density  function  and  the  covariance  functions  will  be 
derived . 

The  characteristic  function  approach  is  to  be  used  in  the 
derivation  of  the  expectation  as  follows: 

E{S)  »  E(S(0  ,cp) }  =  E  ( - i-s  f  [  A(u,v)ej(u0  +  V(p)  dudv 

\  W  J  J 


J A(u,v)  E  |ej(ue+  Vq>)  J  dudv 

=  — f  (  A(u,v)  4>  (u,  v)  dudv 

(2rt)ZJJ  9>V 


(3.32) 


frtfWWMfVs  $r*#****.\  •,.*•***  -  * 
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where 


A(u,  v)  mff  s<e,cp)  e‘J(u0+  de  dcp  (3.33) 

is  the  Fourier  transform  of  S(0  ,<p) ,  and  ^  is  the  characteristic 
function  of  Q  and  <p.  Similarly,  for  higher  order  moments,  the 
expressions  become 


— (//*<»,  v)  eJ(u0+  ^  dud,) 


and 


+  vr»  V1  +  •  •  •  +  vr)  X  duidvi  •  •  •  durdvr  (3.34) 

E  {Sl[e(*l)>  *(*J  s2  [e  {h)'  »(*j)] 

'  7^?  fill  )  **  M  V2-V*2  (ul>vl,u2>v2) 


X  du^  dv^  du2  dv2 


(3.35) 


Thus,  the  moments  and  the  correlation  functions  of  the  scattering 
parameters  can  be  computed,  at  least  in  theory.  From  a  practical 
point  of  view,  however,  this  approach  appears  to  have  the  same  diffi¬ 
culty  as  determining  the  densities  using  the  transformation  of 
variables.  The  digital  computer  makes  the  computations  of  the  moments 
and  the  covariance  functions  possible. 

The  procedure  used  in  the  computational  process  is  basically  to 
reduce  the  integrals  to  summations  using  approximation  techniques. 

The  general  steps  are: 
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(1)  Compute  the  Fourier  integral  A(u,  v)  in  (3.33). 


(2)  Compute  the  characteristic  function  (u,v)  which  for  the 


assumption  of  a  jointly  Gaussian  density  for  the  aspect  angles  is 


given  by 


j(eu  +  cpv^  -  (afu2  +  2p  a. a  uv  +  a2v  j 
*  <u,v)  .  e  '  )e  2  Kb  6  6  *  ‘P  ' 

0|Cp 


(3)  Compute  the  integral  (3.32). 

There  are  practical  considerations  which  influence  the  detailed 
procedures,  which  will  be  discussed  as  they  are  encountered. 

It  was  shown  in  Section  3.2  that  the  conditional  variances  for  0 
and  cp  are  small.  It  is  reasonable,  therefore,  to  consider  only  the 

portion  of  the  0  -  cp  plane  having  a  reasonable  probability  of 

—  —  2  2 
occurring.  Let  0  and  cp  be  the  conditional  means  and  n  and  a  be  the 

-  0  cp 

conditional  variances  of  0  and  cp.  Let  M  be  a  positive  integer.  Then 
the  area  of  the  plane  having  coordinates  such  that  0-Mo.  -  6  -  W 

U 

+  Ho„  and  ra-Mo  ^  cp  ^  cp  +  M  c  will  be  used  for  the  actual 

0  cp  K  ^  v 

computations . 

Now,  let  ST(0,cp)  be  a  periodic  function  which  is  identical  to 
S(0 ,  cp)  in  the  area  defined  above,  i.e., 


-cp-cp-Mo  , 

0  0  cp  K  ^  cp 

and 

ST  ^0  +  2kMoQ  ,  cp  +  2£Ma^  =  ST(0  ,cp) 

for  k,  i  integers. 

The  periodic  function  S^,(0,cp)  instead  of  S(0,cp)  will  be  used  in 
computing  the  expectations.  A  bound  on  the  resulting  error  in  using 
a  finite  M  is  found  in  Appendix  C .  A  reasonable  compromise  seems  to 
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be  M  =  4  if  0  and  cp  are  assumed  to  be  Gaussian.  Since  S^(0,cp)  Is 


periodic,  it  can  be  represented  by  its  Fourier  series. 


00  00 


ST(e,cp)  a  E  E  A(m,n) 


j(nmQe 


(3.36) 


m=  -co  n  =  -oo 


where 


0  Mn 


0  Ma 


Substitution  of  (3.36)  in  the  expectation  results  in  a  discrete  version 


of  (3.32). 


CO  00 


E(S)«E<  E  E  A(m 


m*  -oo  n»  -oo 


00  00  /  \ 

E  E  A(m,n)  <S>  (mu0>  nvQJ  (3.37) 

a  n  22  am  *  r  \  ' 


j  fmu(P  +  nv0a) 
>n)  e  '  ' 


ma  -«  n  =  -oo 


so  the  double  integral  in  (3.32)  becomes  a  double  summation.  The  only 
approximation  made  at  this  point  is  the  substitution  of  S^,(0  , cp)  for 
S (0  ,cp)  as  mentioned  above.  In  practice  only  finite  number  of  terms 
will  be  used  in  the  summation  (3.37).  A  bound  on  the  resulting  error 
is  also  discussed  in  Appendix  C. 

Similarly  (3.33)  and  (3.34)  become,  respectively, 

(  coco  coco 

e  sr  a  E  E  ...  E  E  A  (rai » njJ  •  •  •  A  (®r ,  nr j 

I  m„=-oo  n,=-co  m  =-cc  n  =-«  '  '  '  ' 

11  r  r 


-  v*(“oi,vvou1n  i 


(3.38) 
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Now,  it  is  necessary  to  obtain  the  Fourier  coefficients,  A(m,n). 
These  will  be  estimated  using  the  Fast  Fourier  Transform. 

3.3.2  The  Fast  Fourier  Transform 

The  Fast  Fourier  Transform  (FFT)  is  a  digital  computer 
algorithm  for  fast  and  efficient  computation  of  the  Discrete  Fourier 
Transform  (DFT)[61,  62] .  The  DFT  is  defined  by 

Ar  "  5  xn  exp  (3'40) 

where  A_  is  the  r-th  coefficient  of  the  DFT  and  x  is  the  n-th  sample 
r  n 

of  the  time  series  which  consists  of  N  samples.  The  inverse 
transform  is  given  by 

xn  "  X  *x  eXP(^£S)  •  (3‘41) 

One  of  the  consequences  of  using  the  DFT  is  that  both  the  time 
series  and  the  coefficients  of  the  harmonic  frequencies  can  be  defined 
for  all  integers  (positive  and  negative).  This  results  in  the  familiar 
periodic  form 


A  ®  A  ...  =  A  5 
r  r+N  r+2N 

n  nfN  n+2N 


The  approximation  of  the  Fourier  Transform  or  Fourier  Series  by 
by  the  DFT  is  not  without  its  pitfalls  and  hazards.  A  distortion  due 


to  aliasing  occurs  in  the  frequency  domain  for  real  functions  which 
are  not  sampled  at  sufficiently  high  rates  [62,  63]  . 

In  addition  to  the  aliasing  error,  there  is  an  error  due  to  the 
fact  that  the  FFT  acts  like  a  bank  of  filters,  each  with  a  sin  Nx/sin  x 
response,  where  N  is  the  number  of  samples  used  (Appendix  C) .  Let  Br 
be  the  r-th  Fourier  series  coefficient  of  a  periodic  time  function, 
xT(t)  »  x(t) ,  0  £  t  -  T,  then  the  FFT  estimate,  Ar>  of  Br  is  given  by 

1  2  „  "  r^(  N  )*  sin(m  -  r)ir 


A  a  —  Z  Be 
r  N  m 

m  =  -co 


sin(m  -  r)rt 
sin(m  -  r)^ 


Z  B. 


(3.42) 


The  error  in  estimating  the  Fourier  transform  of  a  non-periodic 
function  due  to  the  "sidelobes"  of  the  filter  is  called  leakage  [64] . 
The  replacement  of  x(t)  by  its  periodic  version  x^,(t),  as  mentioned  in 
the  previous  section,  reduces  this  error  to  the  error  introduced  by  the 
harmonic  components  B^  outside  the  unambiguous  range  of  A^.  This  error 
is  due  to  failure  to  satisfy  the  Nyquist  sampling  rate  for  the  function 
being  analyzed.  In  the  work  here,  the  sampling  rate  is  sufficient  to 
make  this  error  insignificant. 

The  DFT  and  FFT,  can  be  generalized  to  two  dimensions.  The  two¬ 


dimensional  DFT  is  defined  by 

.  N-l  N-l 

A  s*  -4r  Z  Z  X 
r,s  .,2  n  n  n 

*  N  n»aO  nsaO 


exp/^2il£] 

i,n  P\  N 


rm  +  sn) 


(3.43) 


where  A^s  is  the  coefficient  corresponding  the  r-th  harmonic  in  one 
dimension  and  the  s-th  harmonic  in  the  other  dimension,  and  xmn  is  the 
m-th  by  n-th  sample  of  x.  In  this  case  there  are  N  samples  in  each 


)*-r 
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dimension  of  x  and  A.  Just  as  in  the  one  dimensional  case,  A  and 

rs 

xmn  can  defined  so  that  they  are  periodic  functions,  but  in  each 
dimension,  i.e., 


A  ~  A  _ _  —  A  ss  A  ..  ...... 

r,s  r,sfN  r+N,s  r+N.&fN 


Xm,n  “  Xm,nhN  “  xnH-N,n  “  XmfN,nfN 


The  inverse  transform  is  defined  by 

N-l  N-l 

x  =  Z  Z  A 
m,n  r=0  s=0  r,s 


^  2rt  j  (rm  +  sn)^ 


(3.44) 


The  two-dimensional  FFT  is  used  to  obtain  the  coefficients  necessary 
to  compute  the  moments  and  covariances  in  (3.37),  (3.38),  and  (3.39). 
The  characteristic  function  is  computed  directly  by  assuming  that  6 
and  cp  are  jointly  Gaussian  random  processes. 

3.3.3  Practical  Limitations 

A  slight  regression  is  in  order  at  this  point  to  aid  in 
visualization  of  the  problems  encountered  in  actually  implementing 
these  equations.  Equations  (3.37),  (3.38),  and  (3.39)  represent  a 
rather  significant  problem  in  terms  of  the  capacity  of  a  computer, 
i.e.,  in  terms  of  memory  capacity  and  number  of  operations  required. 

By  proper  choice  of  the  origin,  the  characteristic  function  can  be 
made  a  real  function.  The  coefficients  of  the  two  dimensional  trans¬ 
form,  however,  are  generally  complex  functions.  If  N  =»  128,  then 
there  are  16,384  real  numbers  required  for  the  characteristic  function 
and  equally  as  many  complex  numbers  required  for  the  A  in  (3.37).  In 
terms  of  memory  requirements,  the  characterisitic  function  requires 


16,384^q  »  40,000g  words  in  memory  and  the  A  requires  32,768^^  = 


nwwK’fcisssswiws 
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100,000g.  This  represents  a  significant  portion  of  the  resources  of 
even  a  large  scientific  computer.  Even  more  significant,  however,  is 
the  number  of  operations,  or  iterations  implied  by  the  multiple  sums. 

Assume  that  (3.37),  (3.38),  and  (3.39)  are  to  be  implemented  as 
written,  which  is  to  say  by  brute  force.  Equation  (3.37)  would 
require  16,384  iterations,  (3.39)  would  require  268,435,456,  and  (3.38) 
would  require  (16,384)  where  r  is  the  moment  desired.  If  each 
iteration  required  1  microsecond  (an  optimistic  estimate),  a  second 
moment  would  require  about  268  seconds,  clearly  an  unreasonable 
length  of  time  for  one  number.  It  is  possible,  however,  to 
significantly  reduce  this  time,  but  higher  order  moments  are  clearly 
beyond  consideration  using  this  approach. 

It  is  shown  in  Appendix  C,  that  for  the  data  used  here,  the  size 
of  the  arrays  required  for  the  Fourier  coefficients  can  be  reduced  by 
a  factor  of  approximately  4.  This  means  that  the  number  of  iterations 
required  for  (3.37)  is  also  reduced  by  a  factor  of  about  4.  The 
reductions  are  more  significant  for  (3.39).  The  four-dimensional 
characteristic  function  requires  only  about  1/16  of  the  original  memory 
requirements  and  only  about  1/16  as  many  iterations  are  required. 
Further  reductions  in  computation  time  are  possible  due  to  the  rate  at 
which  the  characteristic  function  converges  toward  zero. 

The  third  and  fourth  central  moments  are  also  desired  to 

characterize  the  probability  density  functions.  These  are  more 

appropriately  calculated  using  numerical  integration  techniques.  Use 

2 

of  numerical  integration  provides  the  moments  by  using  only  N 
iterations.  For  rectangular  (Euler)  integration,  the  moments  are 
computed  by 
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N-l  N-l 

S  £  Sr(mA0  ,  nAp)  f  ,  nAr>)  A0  Acp 

n*0  i*0  6,Cf) 


(3.45) 


where 


A  6  m 


2_% 

N 


and 


2M  a 


) 


and  where  the  coordinates  have  been  shifted  appropriately.  The  first 
two  moments  computed  using  rectangular  integration  have  been  compared 
with  those  computed  using  (3.36)  and  (3.38)  with  excellent  results. 

3.3.4  The  Computation  of  the  Statistics 

The  previous  sections  of  this  chapter  have  developed  the 
techniques  for  the  computation  of  the  statistical  characteristics  when 
the  aspect  angles  are  joir.dy  Gaussian  random  processes  and  the 
scattering  parameters  rre  known.  Th^se  techniques  will  now  be  applied 
to  the  scattering  model  developed  in  Chapter  2. 

The  application  to  the  moments  and  the  covariance  of  the  radar 
croc"  section,  S,  has  already  been  discussed  in  Section  3.3.1.  The 
azimuth  and  elevation  errors  given  in  (2.5),  (2.6),  and  (2.13)  will 
now  be  discussed.  Consider  (2.5),  (2.6),  and  (2,13).  The  azimuth  and 
elevation  errors,  AZER  and  ELER,  are  obtained  by  dividing  the  nonradial 
components  of  power  by  the  radial  component  of  power.  This  means  that 
the  correlation  or  covariance  between  RCS  and  AZER  or  RCS  and  ELER 
will,  in  general,  be  small  since  AZER  and  ELER  are  bipolar  functions 
with  small  means.  In  fact 


S  x  AZER 


=  E 


2i 


cos 


(“i '  °j)} 


(3.46) 
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E  |  S  X  ELERj  .  E  |  ^  £ /s^J  £u  cos (  a,  -  0j)J 


.(3.47) 


Therefore,  more  meaningful  variables  are  the  nonradial  components  of 


normalized  power  instead  of  AZER  and  ELER,  i.e.,  the  numerators  of 
(2.5)  and  (2.6),  namely 


M  M 

Z  Z  /s  s' 

i=l  j=l  1  J 


'(“i '  °j )  - 


S  X  ELER 


(3.48) 


MM  _ 

e.j  =  Z  Z  J  S.S.  f,  cos  (a,  -  a,  S  x  AZER  .  (3.49) 

i=l  2  '  1  J' 

The  statistics  of  these  two  functions  and  S  have  been  computed  for 
a  number  of  different  conditions.  Sample  points  on  a  straight  and 
level  flight  were  selected  and  combinations  of  the  variances  and  time 
constants  of  the  Euler  angles  were  used  to  examine  the  moments  and  time 
constants  of  these  scattering  parameters. 

Figures  3.2,  3.3,  and  3.4  are  computer  plots  of  S,  e^,  and  ££, 
respectively,  for  the  region  of  the  0  -  cp  plane  for  1.5677  -  6  -  1.7317 
and  -0.28241  f  .-p  S  -0.12131.  These  data  represent  a  region  of  +4cr 
from  the  mean  in  each  direction  for  0  =  1.6497,  cp  =  -0.20186,  a  » 

0.02050,  and  a  =  0.02014.  If  the  variances  are  reduced  by  1/4  these 

cp 

data  represent  a  region  of  ±8a.  Table  3.1  contains  the  moments  of  the 
scattering  parameters,  S,  and  £2  for  the  two  different  conditions. 
It  is  clear  that  the  moments  of  these  scattering  parameters  are  highly 
sensitive  to  the  statistics  of  the  aspect  angles. 

Chapter  4  presents  data  from  a  number  of  sample  points  on  the 

flight  path  with  different  sets  of  conditions  on  the  Euler  angles. 

The  results  of  the  analytic  solutions  are  compared  with  the  results 
of  simulation  to  verify  the  techniques  derived  in  this  chapter. 
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S  4.  COMPARISON  OF  ANALYTIC  SOLUTIONS  AND  SIMULATIONS 

t 


4.1  General 

A  straight  and  level  incoming  flight  path  with  a  small  crossing 
angle  was  selected  as  being  typical  of  the  conditions  of  interest. 


Three  sample  points  on  this  flight  path  were  selected  at  which  compu¬ 


tations  of  the  statistics  of  the  scattering  parameters  were  to  be 
calculated.  Figure  4,1  depicts  the  flight  path  and  the  sample  points. 
The  coordinates  of  the  three  sample  points  are  (30,  10,  and  2.5  km), 

(20,  4.238,  and  2.5  km),  and  (10,  -1.524,  and  2.5  km).  In  order  to 
investigate  the  sensitivity  of  the  scattering  parameters  to  the 
aspect  angle  statistics,  three  different  combinations  of  Euler  angle 
statistics  were  applied  at  each  sample  point,  and  the  resulting  aspect 
angle  statistics  were  computed  analytically  and  by  simulation. 

Table  4,1  contains  the  data  on  the  resulting  aspect  angle 
statistics  used  for  each  run.  Figures  4.2  through  4.10  show  the 
scattering  parameters  for  each  of  the  three  sample  points  as  functions 
of  the  aspect  angles.  The  limits  on  6  and  9  are  ±4a  for  the 
largest  values  of  and  at  each  sample  point.  Figures  4.2  through 

4.4  represent  S,  e^,  and  ^  at  sample  point  1;  Figures  4.5  through  4.7 

represent  S,  e, ,  anv’  t-,  at  sample  point  2;  and  Figures  4.8  through  4.10 

represent  S,  e^,  and  €2  at  sample  point  3. 

4.2  Simulation  Results 

A  simulation  program  «s  been  written  which  performs  a  Monte  Carlo 
simulation  of  the  time  history  of  the  EM  scattering.  The  simulation 

t-  users'*-. 
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uses  the  same  assumptions  and  models  that  are  used  in  the  analytic 
solutions,  except  that  the  angular  transformation  is  not  linearized 
about  the  mean  aspect  angles.  The  simulation  is  fairly  general  in 
that  the  flight  path,  velocity,  random  motion  statistics,  and  sampling 
interval  are  input  variables.  The  output  data  are  2048  time  samples, 
taken  at  50-millisecond  intervals,  of  each  of  the  scattering  parameters 
and  their  statistics.  The  aspect  angles  statistics  are  also  available 
as  output  data.  The  simulation  permits  an  investigation  of  many 
characteristics  not  conducive  to  analytic  solution.  For  example, 
by  changing  the  velocity  it  is  possible  to  investigate  the  limits  the 
flight  path  and  flight  characteristics  place  on  the  assumptions  used  in 
Chapter  3  regarding  the  wide-sense  stationarity  of  the  aspect  angles. 

Figures  4.11  through  4.19  are  the  outputs  of  the  simulation  for 
run  6.  Figures  4.11  through  4.13  represent  S,  tp  and  e2;  Figures  4.14 
through  4.16  represent  the  autocovariances  of  S,  ep  and  e2;  and 
Figures  4,17  through  4.19  represent  the  cross-covariances  of  S  and  ep 
S  and  €2,  and  ^  and 
4.3  Data  Comparison 

The  statistics  of  the  scattering  parameters  have  been  computed 
analytically  and  from  simulations  for  each  of  the  runs  defined  in 
Table  4.1.  The  means  and  variances  of  the  aspect  angles  are  given  for 
the  linearized  solutions  and  the  simulations  to  indicate  the  general 
level  of  errors  encountered.  Tables  4.2,  4.3,  and  4.4  compare  the 
statistics  of  the  scattering  parameters  for  the  nine  runs  actually 
made.  Figures  4.20  through  4.37  cu»\pare  the  covariance  functions  for 
runs  1,  6,  and  9. 
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It  is  clear  from  the  data  that  the  statistics  of  the  scattering 
parameters  are  very  sensitive  to  the  statistics  of  the  aspect  angles. 
Some  of  the  differences  between  the  analytic  results  and  simulations 
can  be  explained  by  the  apparently  small  differences  in  the  means  and 
variances  of  the  aspect  angles.  More  important  than  the  differences 
are  the  similarities  between  the  sets  of  data.  It  indicates  that  the 
general  approach  will  permit  the  analytic  determination  of  the  statis¬ 
tics  of  the  scattering  parameters  with  some  degree  of  accuracy. 

The  sensitivity  of  the  scattering  parameters  to  the  statistics  of 
the  aspect  angles  emphasizes  the  need  for  actual  data  on  the  flight 
characteristics  of  aircraft.  This  may  become  even  more  important  when 
higher  frequencies  are  used  since  the  lobing  structure  of  the 
scattering  parameters  increases  with  increasing  frequency. 
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Figure  4.4.  Azimuth  Component  of  Nonradial  Power,  e_, 
at  Sample  Point  1  c 
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Figure  4.5.  The  Radar  Cross  Section,  S,  at  Sample  Point  2 


Figure  4.6.  Elevation  Component  of  Nonradial  Power 
at  Sample  Point  2 


Figure  4.7.  Azimuth  Component  of  Nonradial  Power, 
at  Sample  Point  2 
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Figure  4.14.  Autocovariance  of  S  Versus  Delay  Time 


Figure  4.15.  Autocovariance  of 
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Figure  4.23.  Crosscovariance  of  S  and  Run  1 
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Figure  4.24.  Crosscovariance  of  S  and 
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Figure  4.30.  Crosscovariance  of  S  and  €2,  Run  6 
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Figure  4.31.  Crosscovariance  of  and 
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5.  SUMMARY,  CONCLUSIONS,  AND  APPLICATIONS 
5.1  Sumnary  and  Conclusions 

The  objective  of  this  research  has  been  to  develop  techniques  for 
the  analytic  calculations  of  the  statistical  characteristics  of  the  EM 
scattering  parameters  of  airborne  radar  targets.  In  Chapter  2,  a 
deterministic  model  of  the  EM  scattering  is  developed  and  applied  to 
an  ellipsoidal  model  of  a  target  drone.  In  Chapter  3,  techniques  are 
derived  to  compute  the  conditional  statistical  characteristics  of  the 
scattering  parameters.  These  statistical  characteristics  are 
conditional  on  the  target  flight  path  since  the  mean  aspect  angles  are 
dependent  only  on  the  flight  path  with  respect  to  the  radar.  Chapter 
4  compares  the  analytically  computed  statistical  characteristics  at 
three  sample  points  on  a  flight  path  with  results  of  a  Monte  Carlo 
simulation.  At  each  sample  point  three  sets  of  Euler  angle  statistics 
were  applied  to  investigate  the  sensitivity  of  the  statistics  of  the 
scattering  parameters  to  the  statistics  of  the  aspect  angles.  The 
results  indicate  the  validity  of  the  approximations  used,  at  least  for 
the  selected  target  and  flight  path. 

The  prime  objective  of  developing  techniques  for  analytic  calcu¬ 
lation  of  the  statistical  characteristics  of  the  EM  scattering  param¬ 
eters  has  been  accomplished.  The  ability  to  examine  the  statistical 
relationships  among  the  scattering  parameters  will  be  of  great  value 
in  designing  tracking  filters  and  evaluating  system  performance  against 
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other  than  Ideal  targets.  Furthermore,  the  deterministic  model  of  the 
scattering  parameters  may  aid  in  the  understanding  of  the  scattering 
phenomena  and  their  effects  on  radar  systems. 

The  secondary  objective  of  identifying  areas  where  existing  data 
or  theory  does  not  support  the  general  application  of  the  primary 
objective  has  also  been  accomplished.  The  number  of  problem  areas 
exceeded  expectations.  The  ability  to  model  the  EM  scattering  from 
relatively  simple  geometrical  shapes  is  limited.  Most  aircraft  are 
made  up  of  large  numbers  of  not  so  simple  geometrical  shapes  which 
cannot  in  general  be  accurately  modeled  by  quadratic  surfaces,  circular 
cylinders  or  flat  plates.  It  is  essential  that  key  reentrant  structures 
such  as  the  engine  ducts  be  modeled  since  these  represent  primary 
scattering  centers  in  certain  regions  of  aspect  angles,  in  fact,  for 
many  aircraft  they  are  the  largest  backseat ter ing  elements  in  the  entire 
sector  within  45  degrees  of  nose-on.  The  inclusion  of  blstatic  angles 
and  the  effects  of  polarisation  will  add  additional  complexities  to  an 
already  complex  problem. 

The  almost  total  lack  of  measured  data  on  random  aircraft  flight 
characteristics  is  one  of  the  most  surprising  facts  encountered.  These 
data  are  necessary  to  Improve  the  accuracy  of  the  results. 

Measured  glint  data  are  extremely  limited  and  do  not  exist  on 
targets  of  Interest  for  air  defense  systems.  Therefore,  an  evaluation 
of  the  general  accuracy  of  the  calculated  glint  data  is  limited  to 
empirical  judgment. 


'  *<#*-** 
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5.2  Applications  and  Recommenda tlons 

The  results  of  this  research  have  several  applications  in  radar 
tracking  of  aircraft.  Two  of  these  applications  are: 

(1)  The  statistics  of  the  radar  scattering  characteristics  may 
be  used  in  radar  system  simulations  and  evaluations.  Present  Simula* 
tions  and  analyses  either  ignore  the  target  induced  errors  or  assume 
that  they  are  independent  additive  processes.  It  is  shown  here, 
however,  that  the  target  induced  errors  are  not  independent,  but  that 
they  are  correlated  with  the  RCS  and  each  other.  More  accurate  simula¬ 
tions  may  result  from  more  accurate  modeling  of  the  target  induced 
effects  and  the  correlations  among  them.  In  order  to  further 
simplify  the  application  of  the  results  to  radar  simulations, 
additional  work  should  be  directed  toward  the  development  of  dynamical 
models  of  stochastic  processes  having  the  desired  statistics.  Haddad 

[65]  presents  an  approach  for  dynamical  modeling  of  separable  Markov 
processes,  however,  this  may  not  be  sufficient  for  some  of  the  tar¬ 
gets.  More  complicated  techniques  can  be  based  on  the  works  of  Doob 

[66] ,  Frost  [67],  Faurre  [68],  Kailath  and  Frost  [69],  and  Haddad 
and  Valisalo  [ 70]  . 

(2)  The  error  statistics  may  be  applied  to  the  design  of  optimum 
tracking  systems.  The  usual  approach  to  the  design  and  analysis  of 
tracking  filters  is  to  assume  that  the  measured  errors  are  due  only 

to  additive  white  noise.  Typical  examples  of  such  analyses  include 
the  works  of  Mosca  [71],  Sharenson  [72],  and  Urkowitz  [73].  With 
statistical  models  of  the  target  induced  errors  it  will  be  possible  to 
generalize  these  works  to  more  realistic  conditions,  including  non- 
gaussian,  trajectory-dependent,  correlated  noise. 
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Other  applications  include  the  investigation  of  the  statistical 
relationships  of  the  scattering  parameters  for  different  polarizations 
and  different  frequencies.  Basalov  and  Ostrovityanov  [32]  and  Ewell, 
Alexander,  and  Tomberlin  [74]  have  investigated  the  effects  of 
polarization  on  the  radar  scattering  characteristics  of  a  target  and 
the  effects  on  angle  tracking.  More  accurate  statistical  models  would 
be  useful  in  this  arna  to  evaluate  the  possibility  of  using  polariza¬ 
tion  agility  to  aid  in  reducing  tracking  errors.  The  analytic  model 
developed  in  Chapter  2  can  be  generalized  to  include  the  effects  of 
polarization  by  us'ng  the  works  of  Uslenghi  and  Lee  [55],  Moll  and 
Seecamp  [48],  and  Pierson  and  Clay  [56].  Sims  and  Graf  [34,  35]  have 
investigated  the  use  of  frequency  agility  to  reduce  rms  tracking  errors 
due  to  target  induced  errors.  Finally,  a  statistical  analysis  of  the 
frequency  sensitivity  of  target  scattering,  using  techniques  developed 
here,  may  be  used  to  extend  the  work  in  this  area. 

5.3  Additional  Efforts  and  Future  Plans 

This  research  is  a  portion  of  an  ongoing  task  in  the  area  of 
airborne  radar  target  signatures.  Three  tasks  have  been  planned  to 
help  fill  the  gaps  in  data  and  theory  as  applied  to  this  problem. 

Other  tasks  are  under  development  but  require  coordination  and  concur¬ 
rence  of  other  agencies  for  actual  implementation. 

A  contract  effort  is  in  progress  at  the  University  of  Illinois  to 
develop  solutions  to  the  EM  scattering  theory  for  selected  geometrical 
bodies.  These  will  be  used  first  to  develop  a  more  accurate  model  of 
the  scattering  from  the  BQM-34A  target  drone.  Later  these  solutions 
and  others  can  be  used  to  model  the  radar  scattering  from  other  airborne 
targets  of  interest. 


Two  radar  scattering  measurements  programs  are  In  progress  which 
will  provide  data  to  aid  in  validation  and  improvement  of  the  analytic 
solutions  to  the  scattering  problem.  One  measurement  program  is  at 
RATS CAT,  Holloman  Air  Force  Base,  New  Mexico.  These  measurements 
include  glint,  monostatically  and  bistatically,  in  addition  to  the  more 
conventional  RCS  measurements  on  the  BQM-34A  and  BQM-34F  target  drones. 
The  second  measurements  program  is  being  done  by  General  Dynamics, 

Convair  Aerospace  Division,  Fort  Worth,  Texas  on  the  same  target  drones. 
This  program  will  use  short  pulse  technology  and  sophisticated  data 
processing  to  identify  the  location  and  provide  infor  mation  on  individual 
scattering  centers.  It  is  expected  that  these  two  measurements  programs 
will  aid  greatly  in  the  development  of  more  accurate  analytic  modeling 
techniques . 

The  data  desired  on  aircraft  flight  character is  tics  are  essential 
to  the  statistical  modeling  effort.  However,  any  program  to  make  the 
necessary  measurements  will  require  the  assistance  of  the  Air  Force  or 
the  Navy,  or  both.  Development  of  actual  data  requirements  has  not  been 
initiated  due  to  lack  of  time. 

The  work  reported  here  is  of  interest  to  all  servicec.  Closer 
coordination  is  planned  in  the  future  due  to  the  overlapping  interests 
and  the  expense  involved  in  this  type  of  research. 
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d^  m  sin  0  cos  cp 
d..  .  -  sin  0  sin  cp 
d^  ■  cos  0 
^21  a  cos  ®  cos  cp 

^22  “  cos  6  8*n  ^ 
d22  "  “  sin  0 
d^  ■  -  sin  cp 
d^2  “  cos  cp 

d33  -  0  (A. 3) 

The  incident  wave  has  electric  and  magnetic  components  given  by 

t.  i  +  E  .  ?\ejWt  (A .4) 

inc  y  0  inc  Q  cp  inc  cp )  ' 

and 

fflnc  ■  (H0  lac  \  *  \  inc  %)  <A-5> 

where 

?lnc  -  Re  (finc  ><  O 

The  components  o£  the  electromagnetic  wave  scattered  from  the 
i-th  element,  having  cross  section  S^,  are  expressed  as 

t  ,  »  -  JIT  t. 

ref  v  i  inc 


(A  .6) 
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ttre£--^  fi,= 


(A. 7) 


where 


ref 


•  Si  ?inc 


The  units  of  P  r  are  watts  and  those  of  P.  are  watts/m  . 
ref  inc 

A. 2  Geometrical  Computations 

Let  the  source  and  observation  point  Q  be  located  at  (x,  y,  z)  or 
(R,  0,  9)  and  the  i-th  scatterer  be  located  at  (x^,,  y^,  z^)  with 
respect  to  the  target  centered  coordinate  system.  Then  the  observation 
point  is  at  (x  -  x^,  y  -  y^,  z  -  z^)  or  (R^,  e cp^)  with  respect  to 
the  shifted  coordinate  systems  centered  at  the  i-th  element,  and  has 
unit  vectors  I?ri,  in  the  local  polar  coordinate  system 

centered  at  the  i-th  element.  The  relations  between  the  coordinates 


are  given  by 


Rt  -  [(x  •  *i)2  +  o  -  yi y  +  (*  -  2i)z  j 

(R) 


arctan 


(A. 8) 


(A  .9) 


and 


a  arccos 


-ft 


(R) 


(A.  10) 


The  Cartesian  unit  vectors  i,  j,  k  are  related  to  the  polar  unit 
vectors  u^,  u^,  T?  ^  at  the  observation  point  by  the  transformation 


r-»  1 
u  , 

L-il 

ri 

— > 

j 

U0i 

t 

T?  . 

m  m 

<pi_ 

(A. 11) 


where  is  obtained  by  substituting  0i  and  cp^  for  0  and  cp  in  D. 


The  reflected  E  field  from  the  i-th  element  may  be  written  as 

*i  ■  \i  v +  V  v 

where  the  radial  dependence  is  neglected.  Hence,  the  transformations 
(A. 2)  and  (A. 11)  result  in  the  following  expressions  for  the  components 
of  the  reflected  field  in  the  original  system  are  given  by 


E'  u 
ri  r 


ei  ue 

e'  u 

.  q*  <p. 


[D][DIlJ  Eei  3 


E  .  u  . 
q>i  cpi 


(A. 12) 


If  it  is  assumed  that  R  >>  ,  then  8 1  «  6,  cp^  a*  q>,  and 

Rj,*  R,  then  the  components  of  the  E  field  in  (A.  12)  become  approxi¬ 


mately 


**  Egi  008  0  c°8  9  “  cos  0  sin  <p  -  z^  sin  0^ 
+  ^-Xi  sin  :p  +  y£  cos  cp  ) 


E0i  *  E0i 
%i  *  Ecpi 


(A.  13) 


The  relative  phase  difference  between  the  reflected  electromag¬ 
netic  waves  of  two  scatterers  is  given  by 


°-TiE 


(A.  14) 


where  AR  is  the  difference  in  range.  The  phase  shift  for  each 
scatterer  must  be  computed  so  that  the  vector  fields  can  be  properly 
summed.  If  the  target  center  is  used  as  the  reference  point,  then  for 
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the  i-th  scatterer  AR^«  R  -  R^  «  -  sin  0  cos  «p  -  y^  sin  0  sin  q) 
-  z^  cos  0  so  that  the  fields  observed  at  point  Q  from  the  i-th 
element  have  a  time  dependence  of  e^  +  . 

A. 3  Total  Field  Summation 
Let 


it  ■  P  uy  ■  Re  (1?  x  if*)  if. 

be  the  incident  power  on  the  target  array  where 

1?  »  E  (cos  5  tL  +  sin  5  u  | 

\  0  <j>/ 

1?  -  H  (-  sin  6  u„  +  cos  6  I?  ^ 

V  0  <P/ 


(A.  15) 


and  6  is  the  polarization  angle  of  the  E  field  referenced  to  the  -u, 

d 

vector.  Then  the  electromagnetic  fields  reflected  from  the  i-th 
element  and  observed  at  Q  are  given  from  (A. 13)  by 


Erec  t  ’  I  r*l  [(*P  ‘u  -  ^  hi)\ 

+  cos  8  tf,  +  sin  6  if  I 

0  cp  J 

..  H  r-[/  sin  8  r  .  cos  8  ,  \ -» 

Hrec  i  R  v  Si[\“  “T"  fli  +  ~R~  f2i )  ur 


(A. 16) 


-  sin  8  u  +  cos  8 

0 


(A. 17) 


where 


f  ^  ■  “  xj_  cos  ®  cos  V  “  cos  0  sin  q>  *  sin  0 


and 


-  sin  cp  +  y^  cos  cp 
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Hence,  the  total  fields  at  Q  due  to  the  entire  target  array  are  given 
by: 


■W  -I  [(*H  *11  +  *2l)  +  co.  8 

1  j(wt  +  a,  +  *) 

+  sin  8  u  e  v  1  7 

<pj 

-roc  *  f  |1/Ti[(-JT£  *11^*21)“.- 

j(wt'+  a£) 


e 


(A.  18) 


cos  8  u 


6 


+  cos  8  t?  e 


t  ' 

<P 


(A.  19) 


The  components  of  power  received  at  Q  are  found  from  (A.l).  Thus, 
the  iir  component  of  received  power  is  given  by 


Er  rec  “  Ee  (Ee  rec  Hcp  rec  “  Ecp  rec  rec) 

p  M  M  __ 

■  —s  £  £  /Ts,  cos  fa,  -  a.  j  .  (A 

RZ  i-1  j-1  "  J  \  1  J  / 


.20) 


Similarly  the  u^  and  components  of  received  power  are  expressed, 
respectively,  by 


cp  rec 


Re  (e  H*  -  Ert  H*  \ 
\  r  rec  e  rec  e  rec  r  rec  / 

p  M  M  _ 

-o  £  £  y  s  s  f  cos  fa.  -  a.  ) 

k  i-l  j-1  1  j  \  1  J  / 


(A. 21) 


and 


Pn  -  Re 
6  rec 


( E  H*  -  E  H*  ^ 

^  q>  rec  r  rec  r  rec  <p  rec  J 

p  M  M  _ 

—z  £  £  y  s.s.  f . .  cos  (a.  -  a,  j 

RJ  i-l  j-l  1  J  11  '  1  J  / 


(A.  22) 
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A. 4  The  Scattering  Parameters 


The  RCS,  S,  Is  defined  by 


(A.  23) 


with  both  Finc  and  Pre^  being  referred  to  the  target,  that  is  with  the 

2 

radial  dependence  of  power  removed.  has  units  of  watts/m  and 

P  .  has  units  of  watts.  In  this  case, 
ref 


M  M 


S  -  —  p 

P  r  rec 


E  E  J  s's,’  cos  (a.  -  a.  ^  .  (A . 
i-1  j-1  13  V  1  3  / 


The  u^  error,  in  linear  units  at  the  target,  is  given  by 


R  P  .MM 

cd  rec  1  v  v 


cp  P. 


r  rec 


mlS  >1  7’5i*I  £2i  c°s  (“l  -  aj  ) 


(A  .25) 


The  u  error,  e  ,  measured  in  linear  units  at  the  target,  is  given  by 

U  U 

e°  ’  ”*777  "  ®  i^i  fli  cos  (“i  ’  °j)  - 


(A.  26) 


The  target  phase,  referenced  to  the  target  center,  is  given  by 


PHAS  ■  arctan 


EyiT  sin  a. 

i-1  1 _ _ 

M 

E  J  s.  cos  a. 
i-1  1  x 


.  (A.  27) 


These  equations  are  used  in  Chapter  2  in  the  analytic  modeling  of 


the  scattering  parameters  from  a  target. 
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APPENDIX  B.  STATISTICS  OF  THE  ASPECT  ANGLES 


Let  (Rq,  6q,  vPq)  be  the  polar  coordinates  of  the  target  In  the 
ground  fixed  coordinate  system,  and  (R,  0,  cp)  be  the  coordinates  of 
the  origin  of  the  ground  fixed  coordinate  system  in  the  target  fixed 
coordinate  system.  Since  Rq  «*  R,  only  the  aspect  angles  need  to  be 
considered  in  the  analysis.  From  (2.1)  and  (3.3) 


where  q!  ■  (p  -  a.  The  angles  9q  and  cp^  are  dependent  on  the  flight  path 
only  and  therefore  they  can  be  considered  as  deterministic  for  the 
calculation  of  the  conditional  joint  density  of  6  and  cp.  Therefore, 
(B.l)  and  (B.2)  may  be  considered  as  functions  of  the  random  variables 
a,  p,  and  y ,  and  p  and  y  respectively.  It  is  possible  to  solve  (B.2) 
for  p  and  y  as  functions  of  Q  and  $  as  follows. 


Let 


then 


Let 


then 


r0  cos  0O  “  8in  e0  cos  Vo 

rQ  sin  pQ  -  cos  eQ  ,  (B.3) 


sin  9  cos  <p  ■  -  Tq  cos  Pq  cos  p  -  sin  Pq  sin  p 

■  “  rQ  COS  (p  -  P0)  .  (B.4) 


/  ,  2  „  7~2~Z~  2~ 

r  sin  6  sin  cp  +  cos  6 
r  cos  7q  »  sin  9  sin  <p 

r  sin  7q  »  cos  6  ,  (B.5) 


sin  Qq  sin  cpQ  *  -  r  cos  y^  cos  y  +  r  sin  y q  sin  y 

-  -  r  cos  (7  +  7q)  .  (B,6) 

Therefore,  from  (B.4)  and  (B.6), 


cos  (p  -  p0; 

v  sin  6 

cos 

«v 

(B.7) 

cos  (7  +  70) 

sin  9q 

r 

sin 

Vo 

(B.8) 

tan  pQ  - 

cos  0Q 

(B.9) 

sin  6q  cos 

~Q 

tan  7q  « 

cos  d 

(B.10) 

sin  6  sin 

“W 

<P 

• 

Hence,  p  and  y  can  be  written  explicitly  as 
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then 


Let 


then 


rQ  cos  0O  ■  sin  0Q  cos  cpQ 

rQ  sin  pQ  -  cos  0Q  ,  (B.3) 

<V 

sin  0  cos  q>  ■  -  Tq  cos  p^  cos  p  -  Tq  sin  p^  sin  p 

-  -  rQ  cos  (p  -  pQ)  .  (B.4) 

y2  2  esj  2 

sin  0  sin  ip  +  cos  0 
r  cos  ■  sin  0  sin  <p 

r  sin  7^  »  cos  0  ,  (B.5) 


sin  0q  sin  cpQ  ■  -  r  cos  y^  cos  y  +  r  sin  y q  sin  y 

-  -  r  cos  (y  +  .  (B,6) 


Therefore,  from  (B.4)  and  (B.6), 


0s  -  »0 

\  sin  0 

’  ro 

cos 

<v 

JG 

(B.7) 

(r  +  70) 

sin  0Q 

r 

sin 

fo 

(B.8) 

cos  en 

t«n  ^0  - 

sin  0 q  cos 

* 

(B.9) 

tan  7q  . 

cos  0 
sin  0  sin 

TV 

<P 

• 

(B.10) 

Hence,  p  and  y  can  be  written  explicitly  as 
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p  «  arccos 


y  m  arccos 


sin  0  cos  %  ^  t 
/  sin  0Q  sin  cp0  \ 


+  arc tan 


arctanl 


/  cos  eQ  v  ^ 
\8in  a0  cos  cpg  J  (B* 

m(rt«T.ln  ?)  <BJ 


where  r^  and  r  are  defined  in  (B.3)  and  (B.5). 

The  joint  probability  density  function  of  0  and  9  is  found  by 


the  transformation 


f„(M 


where  p  and  y  are  expressed  as  in  (B.ll)  and  (B.12),  and  J  is  the  Jaco¬ 
bian  of  the  transformation.  Since  a  is  an  additive  factor  to  9,  the 
joint  probability  density  function  for  0  and  9  is  found  by  convolving 
the  joint  density  function  for  0  and  9  with  the  density  function  for  a 
in  the  9  variable  only. 

The  magnitude  of  the  Jacobian  of  the  transformation  is  given  by 

I  T 1  d0  dq>  & 

where  the  partial  derivatives  are  found  by  implicit  differentiation  to 


cos  y  cos  9 


^0  ,  ~ 
^  =  sin  9 


Qjj. 

-  sin  7  -  cot  0  cos  7  sin  q5 


cot  0  cos 
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resulting  in  the  expression 

/  1  -  »1»2  e0  sin2  fg 

|j|  *y  - - - -  ■ cos  <p  • 

sin  6 

The  transformation  from  p-y  space  to  fi-'cp  space  is  one-to-one  only 
if  p  and  y  are  restricted  to  the  region  0  *  p  ^  *  ,  and  0  S  y  <  2«. 

The  quadrant  for  p^  and  y^  can  be  found  by  considering  the  signs  of 
the  numerator  and  denominator  in  (B.9)  and  (B.10).  The  quadrant  for 
( y  +  y q)  and  (p  -  p^)  are  found  by  considering  (B.17)  as  well  as  (B.8) 
and  (B.ll). 

Figures  B.l  and  B.2  depict  the  joint  probability  density 
functions  for  0  and  and  0  and  q>  for  6  ■  1.6497  and  q>  ■  -0.20186, 
i.e.,  at  sample  point  1  as  defined  in  Chapter  4.  Table  B.l  compares 
the  moments  calculated  by  numerical  integration  of  the  joint  probabil¬ 
ity  density  function  for  6  and  cp,  and  the  moments  calculated  from  the 
linear  approximations  of  (3.19)  and  (3.20). 
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APPENDIX  C.  ERROR  ANALYSIS 


C.l  General 

It  was  pointed  out  In  Chapter  3  that  the  process  described  lor 
computing  the  statistical  characteristics  introduces  errors  due  to  the 
assumptions  made.  In  the  computation  of  the  four  moments  by 
numerical  integration  there  are  errors  introduced  at  two  points.  In 
the  computation  of  the  covariance  functions  using  the  characteristic 
function  method,  there  are  errors  introduced  at  three  points.  This 
Appendix  estimates  these  errors  or  their  level  of  significance. 

C.2  Computation  of  the  Moments  Using  Numerical  Integration 

There  are  two  sources  of  errors  in  the  computation  of  the  moments 
by  numerical  integration  as  described  in  Chapter  3.  One  error  is  due 
to  truncation  of  the  data  and  the  probability  denisty  function.  The 
other  error  is  due  to  the  use  of  the  rectangular  (Euler)  technique  for 
numerical  integration.  Let  x  and  y  be  normal  random  variables  with 
joint  probability  density  function  f(x,y)  and  let  g(x,y)  be  a  function 
whose  moments  are  to  be  calculated.  Both  functions  are  truncated  at 


-4 a  and  +4c  from  the  means  in  each  dimension  and  then  sampled  at  N 
points  in  each  dimension.  The  moments  of  g(x,y)  are  computed  by 

I  «  N  r 

E  g  (x,y) >  ■  E  Z  g  (nAx,mAy)  f(nAx,mAy)  AxAy 
)  i*l  n*l 

(C.l) 

where  the  origin  is  shifted  to  the  (-4a,  -4a)  point  and  Ax  and  Ay 
are  the  sample  intervals.  If  N  is  sufficiently  large,  and  g(x,y)  is  a 
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relatively  smooth  function,  then  the  accuracy  obtained  by  (C.l)  is 
adequate.  The  error  is  computing  the  expectation  of  a  constant 
function  over  the  region  described  using  128  data  points  in  each 
dimension  is  about  0.013  percent. 

The  first  two  moments  computed  this  way  have  been  compared  to  the 
results  using  (3.37)  and  (3.39)  with  excellent  results.  As  was 
demonstrated  in  Chapter  4,  excellent  agreement  with  the  results  of 
simulations  was  also  obtained. 

C.3  Computation  of  the  Covariance  Functions 

The  computation  of  the  covariance  functions  using  a  discretized 
version  of  the  characteristic  function  method  introduces  errors  at 
three  points.  The  one-dimensional  case  will  be  considered  first.  The 
errors  in  computing  the  expectation  of  a  function  g(x)  are  introduced 
due  to  the  following  assumptions: 

(1)  The  replacement  of  the  actual  function  g(x)  by  a  periodic 
function  g^(x). 

(2)  The  use  of  the  Fast  Fourier  Transform  (FFT)  to  estimate  the 
Fourier  series  coefficients  for  the  periodic  function  g^(x). 

(3)  The  use  of  a  finite  number  of  the  Fourier  series  coefficients 
in  the  summation  of  the  expectation  of  gT(x) . 

The  actual  function,  g(x),  is  replaced  by  a  periodic  function, 
gT(x),  such  that 

gT(x)  -  g(x) 


for 
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gj(x)  «  g^x  +  kT) 


where  the  period  T  is  defined  to  be  2M<j  where  a  Is  the  standard  devia¬ 
tion  of  the  probability  density  function  f(x)  which  is  assumed  to  have 
zero  mean. 

The  first  error  in  computing  the  expectation  (3.39)  is  introduced 
by  the  approximation 


E  jg(x)}  s  E  {gT(x)} 
This  error  may  be  bounded  as 


(C.2) 


(  \ 

2 

El" 

E 

|g(x)  -  gjCx) j 

s  / 

-00 

g(x)  -  gT(x) 

/ 


T 

2 


g(x)  -  gT(x) 


f(x)  dx  S  4  sup 


g(x) 


erfc 


f(x)  dx 

(4)  • 

(C.3) 


For 


so  that 


^  ■  4o  ,  erfc  »  3  x  10 


E1  £  1.2  x  10“4  sup  |  g(x) 


(C.4) 


The  periodic  function  g^(x)  may  be  represented  by  its  Fourier  series 

(j  <c.5) 

so  that  its  expectation  is  given  by 


gT(x)  -  Z  Am  exp 
m  «  -co 


E 


M 

Z  A  exp 
m 

m«  -oo 


(C.6) 


where  4>x(v)  is  the  characteristic  function  of  the  random  variable  x. 
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‘  The  second  error  encountered  In  computing  the  expectation  is  in 
estimating  the  Fourier  series  coefficients  Ar  by  using  the  FFT.  From 
(3.40),  the  estimate  of  the  n-th  coefficient  is  given  by 

f  1  ..  .  .  /  .  2rt  kn  Ax\ 

AQ  -  jj  8T(kAx)  exp  1 - )  (C.7) 

which  with  (C  .5)  results  in 


.  N-l  eo 

I  2  2 

k»0  n*~  » 
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2  rt  km  Ax^ 
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kn  Ax  i 
T  ) 


-  z 

co 
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The  resulting  error  is  due  to  aliasing  caused  by  failure  to  meet 
the  Nyquist  sampling  criterion.  Furthermore,  it  only  yields  a  finite 
number  of  coefficients. 

The  third  source  of  error  is  in  taking  only  a  finite  number  of 
terms  in  the  summation  (C.6)  namely, 

E  jgT(x)}  s  Z  A  ®  (n)  (C.9) 

1  1  '  n—  k  n 


where 


%  <■>)  - \ (t2)  -i(f) 


(C.10) 
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A  bound  on  the  resulting  error  may  be  obtained  as 


E, 


‘  A  ®  (n) 
|n|>k  n  x 


S  2  sup 
n  >k 


£  2  sup 

A 

n 

n  >k 

U 

E  $  (n) 

nxkfl 


exp 

■ 

H^)2] 

1  -  exp  ^ 

«  (k+  1)\ 

M  / 

.  (C.ll) 


For  M  ■  4  and  k«  8  the  resulting  bound  is  equal  to 


,-H 


E,  £  3.10  sup  |A 
n  >k 


(C.12) 


The  above  analysis  for  the  third  error  source  does  not  apply  to  a 
joint  characteristic  function  of  strongly  correlated  random  variables. 
The  normalised  two-dimensional  characteristic  function  for  two  jointly 
Gaussian,  random  variables  is  given  by 


2  2 
+  2ruv  +  v 


)] 


$  (u,v)  »  exp  |  (u^ 
where  r  is  the  correlation  coefficient  of  x  and  y.  When  r  ■  1, 

IV 

G>  (u,-u)  ■  1 


so  that  the  argument  for  limiting  the  sum  as  in  the  one-dimensional 

case  does  not  hold.  However,  when  |U+  rv|>8,  the  characteristic 

.9 

function  is  less  than  2  x  10  .  Therefore,  only  about  8  terms  on  each 

side  of  the  line  u  *»  -  v  sgn(r)  need  be  taken  in  the  summations. 

Now  consider  the  spectrum  of  the  transform  of  the  scattering 
parameters  under  study  here.  It  can  be  shown  that  99  percent  of  the 
power  is  contained  in  about  half  of  the  total  spectrum.  Therefore, 
it  is  possible  to  limit  the  summations  to  a  narrow  band  about  the 


line  defined  above.  In  practice,  it  is  possible  to  obtain  excellent 
results  using  even  fewer  components  than  those  containing  99  percent 
of  the  power.  Similarly  the  aliasing  error  represented  by  (C.8)  can 
be  neglected. 


